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UNIVERSITÄT HANNOVER 
INSTITUT FÜR MECHANIK 


DIPLOMARBEIT 








für 


Herrn cand. mach. Dieter Lindemann, Mätr, Ne. 732 297 


In dieser Diplomarbeit soll die Nachrechnung der räumlichen Poten- 
tialströmung in Rotoren von Windkraftanlagen mit horizontalen 
Drehachsen durchgeführt werden. Hierzu sind die folgenden Punkte 
zu bearbeiten: 


1. Eine Einarbeitung in die Literatur über Auslegungsmethoden für 
Windkraftanlagen ist durchzuführen. Daneben sollen auch Geome- 
trie- und Referenzdaten ausfindig gemacht werden, um die eige- 
nen Ergebnisse werten und vergleichen zu können. 


2. Zu einem bestehenden FORTRAN-Rechenprogramm (Panelverfahren) 
ist ein Modul zur Bereitstellung der Rotorblattgeometrie zu 
entwickeln, das flexibel gestaltet sein soll, um mit wenig 
Aufwand eine breite Klasse von Windkraftrotoren untersuchen zu 
können. Zur Erfüllung der Kutta-Joukowski Bedingung an der 
Profilhinterkante ist ein Nachlaufmodell zu erarbeiten, das 
den Gegebenheiten einer realen Rotordurchströmung möglichst 
nahe kommen soll. 


3. Für die Umströmung des Rotors bei achsparalleler Anströmung 
(stationäre Strömung im Relativsystem) sind Parameterstudien 
durchzuführen und zu diskutieren. Dieser Vergleich von theore- 
tischen Kennzahlen mit den in der Literatur üblicherweise ver- 
wendeten Größen soll die Einsatzbreite der verwendeten Theorie 
demonstrieren. 


4. Schließlich ist die Berechnung von instationären Strömungsfel- 
dern theoretisch vorzubereiten, um z.B. Rotoren mit geneigten 
Drehachsen (Wagner-Rotor) oder die Rotor-Turm Wechselwirkung 
berechnen zu können. 


Die Ergebnisse der Arbeit sind in einem Bericht darzustellen. Die 
neu entwickelten Programmteile sind anwenderfreundlich zu gestal- 


ten. 


Der Abgabetermin der Arbeit ist der 19. Februar 1988. 


(Prof. Dr.-Ing. I. Teipel) 
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oo Punkt im Unendlichen 


1 Einleitung 


Die Geschichte der Entstehung des Lebens vom Urknall bis hin zum Menschen kann als 
eine Geschichte der Suche nach Energiequellen gesehen werden [1]. Auch die Geschichte der 
Menschheit selbst kann unter diesem Gesichtspunkt betrachtet werden. Die Windenergie 


spielte dabei eine größere Rolle, als es uns heute erscheinen mag /?.. 


Neben der Nutzung der Windenergie zum Vortrieb, wie auf Segelschiffen. wurde die Wind- 
energie an Land von Windmühlen vor allem zum Mahlen von Getreide und zum Wasser- 
pumpen eingesetzt. Dabei hatte sich eine Bauart unter den Windmühlen durchgesetzt: 
die Holländer Windmühle. Die Holländer Windmühle hat ein Windrad mit horizontaler 
Achse. Der Rotor nutzt die Energie des Windes über die am Flügel erzeugte Auftriebskraft 
aus. (Auftriebsprinzip). Auch heute ist es prinzipiell noch genau diese Bauart, die unter 
den Windkraftanlagen am stärksten verbreitet ist. Aus diesem Grunde sollen im Rahmen 
dieser Arbeit auch nur Windkraftanlagen mit Horizontalachsrotor behandelt werden, die 


nach dem Auftriebsprinzip arbeiten. 


Andere Bauarten von Windkraftanlagen sind in-[3] und [4] beschrieben. Die Stellung der 
hier behandelten Horizontalachsrotoren unter den Windenergiekonvertern kann man Bild 
1.1 entnehmen. Die Horizontalachsrotoren teilt man in erster Linie nach ihrer Schnel- 
laufzahl ein. Die Schnellaufzahl ist das Verhältnis aus Umfangsgeschwindigkeit an der 
Flügelspitze zur Anströmgeschwindigkeit. Langsamläufer mit niedriger Schnellaufzahl er- 
kennt man an einer großen Flügelzahl und einer großen Flügeldichte. Schnelläufer hingegen 


haben eine kleine Flügelzahl. 


2 Aerodynamische Berechnung von Horizontalachsrotoren 


2.1 Geschichtliches zur aerodynamischen Berechnung 


Der Bau von Windmühlen wurde früher als Handwerk betrieben. Erfahrungen wurden 
über Generationen weitergegeben. Dabei gelang es im Laufe der Zeit, die Windmühlen in 
einer gewissen Weise zu optimieren. Damals kam es mehr darauf an, zu wissen, wie eine 


Windmühle gebaut werden muß, als zu wissen warum sie funktioniert. 


Doch schon zu Beginn des 18. Jahrhunderts wurde die Windenergie Gegenstand theoreti- 
scher Untersuchungen. Namhafte Wissenschaftler wie Pitot, Bernoulli, d’Alembert, Euler 
und Coulomb beschäftigten sich mit der Frage der Windenergie [5]. Ihre Erkenntnisse 
hatten jedoch kaum Einfluß auf die Konstruktion der Windmühlen. Ein Grund dafür mag 
wohl gewesen sein, daß die Forscher und Praktiker nicht zusammen arbeiteten. Coulomb 
hatte z.B. das Problem, daß ihm nicht erlaubt wurde, mit einer der 50 Windmühlen in 


seiner Umgebung Versuche durchzuführen. 


Bernoulli machte im Jahre 1738 darauf aufmerksam, daß die relative Geschwindigkeit des 
Windes am Windmühlenflügel berücksichtigt werden müsse. 1742 machte sich Maclaurin 
Gedanken zum optimalen Anstellwinkel der Flügel. Er erkannte, daß sich die Strömungs- 
bedingungen über dem Radius des Flügels ändern und betrachtete daher einzelne Ele- 
mente des Flügels. 1864 faßte Neumann die gewonnenen Erkenntnisse in seinem Buch 


”Die Windmühlen - Ihr Bau und ihre Berechnung” [5] zusammen. 


1891 begann der dänische Staat, eine Versuchsmühle zu finanzieren. Die Mühle stand unter 
der Leitung von La Cour. Die Ergebnisse seiner Forschung wurden in Deutschland im 
Jahre 1905 veröffentlicht [6]. La Cour erkannte u.a., daß bei der Berechnung der Leistung 
einer Windkraftanlage nicht von der Flügelfläche auszugehen ist, sondern daß vielmehr die 
von den Flügeln bestrichene Kreisfläche von Bedeutung ist. Er stellte fest, daß bei einer 
bestimmten Windgeschwindigkeit die Drehzahl eines Windrades mit steigender Flügelzahl 


und steigender Rotorblattflächendichte abnimmt. 


Zu Beginn des 20. Jahrhunderts wurde mit zunehmenden Erkenntnissen auf dem Gebiet 
der Strömungsmechanik auch die Windenergie immer mehr zum Gegenstand internationa- 
ler Forschung. Es wurde versucht, die auf anderen Gebieten schon vorhanden Erkenntnisse 
auf die Windenergie anzuwenden. Pfleiderer übertrug die Berechnungsmethoden von Hy- 
draulischen Strömungsmaschinen auf die Windräder [7]. Erfolgreicher waren jedoch die 
Ansätze, die von den Erfahrungen der Propellerberechnung ausgingen, wie dies bei Betz 
[8} und Glauert [9] der Fallist. Man kann heute wohl feststellen, daß die moderne Windrad- 
berechnung auf diese beiden Personen zurückgeht. Dabei sind im deutschen Sprachraum 


die Betzschen Ansätze verbreitet, während man im englischen Sprachraum von Glauerts 


Rechnungen ausgeht. 


1926 erschien das Buch ”Wind-Energie” von Albert Betz [8]. Betz definierte Wirkungs- 
grade und Kennzahlen, die auch heute noch in dieser Form gebräuchlich sind. Er errechnete 
wohl als erster den theoretisch maximalen Anteil der Leistung, der einer Luftströmung ent- 
zogen werden kann: 16/27 = 59.3%. Diese Zahl wird in der Literatur auch als Betzscher 
Wirkungsgrad geführt. Die maximale Leistung kann dem Wind theoretisch dann entzogen 
werden, wenn die Windgeschwindigkeit um 2/3 auf 1/3 abgebremst wird. Würde die Luft 
stärker abgebremst, sagen wir bis zum Stillstand, dann könnte kein Massenstrom mehr das 
Windrad durchströmen. Den höchsten Massenstrom erhielte man, wenn die Luft mit ih- 
rer ungebremsten Ausgangsgeschwindigkeit das Windrad durchströmen würde, dann wäre 
jedoch aufgrund der fehlenden Abbremsung wiederum keine Leistungsausbeute möglich. 
Betz beschreibt weiterhin, wie auf iterativem Wege aus der Geometrie des Windrades 
und den aerodynamischen Beiwerten der Flügel die Abminderung der durch das Windrad 
strömenden Luft berechnet werden kann ([8] 5.37). Bei der Herleitung seiner Gleichungen 
vernachlässigt Betz jedoch den Drall der abfließenden Strömung und damit die im Nachlauf 


des Windrades verlorene Rotationsenergie der Strömung. 


Der Drall wurde dann 1935 von Glauert berücksichtigt [9]. In Deutschland war es Ulrich 


Hütter, der die Betzschen Gleichungen unter Berücksichtigung des Dralls in seiner Disser- 


tation von 1942 |10] für einzelne Radelemiente neu herleitete. Seine Ergebnisse wurden in 


"Hütte Band IIA” [11] einem großen Leserkreis zugänglich gemacht. 


F. X. Wortmann entwickelte spezielle Laminarprofile für Windkraftanlagen. Seine Profile 
erkennt man an den Buchstaben ”FX”, mit denen die Profilkennzeichnung beginnt. Der 
Windrad-Profilkatalog von Miley [12] enthält neben anderen Profilen auch Wortmann- 
Profile für Windräder. Wortmann beschäftigte sich weiterhin mit der Gestaltung nach- 
giebiger Nabenkonstruktionen (/4] S.80) und war an der aerodynamischen Gestaltung der 


modernen Großwindanlagen beteiligt. 


Heutzutage werden Windradberechnungen mit Computern durchgeführt, die es erlauben, 
die zur Berechnung nötigen Iterationen schnell durchzuführen. Hütter sowie Wilson und 
Lissaman [13] haben EDV-geeignete Rechenprogramme erstellt, die auf der im folgenden 
beschriebenen Blade-Element-Theory beruhen. Rechenprogramme dieser Art werden von 


der Industrie zur aerodynamischen Rotorauslegung genutzt [14]. 


Die neuesten Unternehmungen versuchen, die Strömung im Windrad mit Computerpro- 
grammen zu berechnen, die auf der Potentialtheorie beruhen. Man bekommt hier funda- 
mentale Probleme, dadurch, daß die Strömung als reibungsfrei angenommen werde muß. 
Als Eingabe benötigen diese Programme jedoch lediglich die geometrischen Daten des 
Windrades. Auftriebs- und Widerstandsbeiwert der Flügelprofile werden nicht benötigt. 


Es ist weiterhin möglich, ein vollständiges Bild des gesamten Strömungsfeldes zu erhalten. 


2.2 Die Blade Element Theory 


Die am häufigsten angewandte Methode zur Berechnung von Horizontalachsrotoren ist 
eine Berechnungsmethode, die vom Verhalten eines einzelnen Radelementes ausgeht. Bild 
2.1 zeigt ein derartiges Radelement eines Horizontalachsrotors. Zur Berechnung wird das 
Windrad in R/Ar konzentrische Radelemente aufgeteilt. In englischen Texten wird die- 
ses Berechnungsverfahren "Blade Element Theory”, ”Vortex Theory” oder auch "Strip 


Theory” genannt. 


Um gleiche Anströmbedingungen an jeder Stelle eines Radelementes zu haben, müssen 


zwei Bedingungen erfüllt sein: 
1. Ar muß hinreichend klein gewählt werden. 


2. Die Windgeschwindigkeit am Flügel darf nicht von der momentanen Stellung des 


Flügels während seiner Drehung um die Windradachse abhängen. 


Oft ist die zweite Bedingung nicht erfüllt: 


1. wenn die Grenzschicht der Luftströmung über der Erdoberfläche berücksichtigt wer- 


den muß, wie bei großen Windrädern. ! 


. wenn aus konstruktiven Gründen bei einem sogenannten Horizontalachsrotor die 


[DD 


Achslage nicht genau horizontal gewählt wurde, sondern um 5°...15° aufgestellt 


wurde, wie dies in der Praxis häufig vorkommt. 


In diesen Fällen empfiehlt es sich, die Radelemente über dem Umfang in 360°/AO® Teile 
aufzuteilen. Wird nun AO hinreichend klein gewählt, so kann man für jedes Radelement 


wieder von konstanten Anströmverhältnissen ausgehen. 


! Einfacher als eine Unterteilung des Windrades über dem Umfang in 360° / AO Teile ist eine in [2] 
angegebene näherungsweise Methode, mit der ein Korrekturfaktor berechnet wird. Dieser Korrekturfaktor 
ist abhängig von der Ausbildung der Grenzschicht und der Geometrie des Windrades. Der Korrekturfaktor 


ist in [2] in Tabellen gegeben. 





2.3 


Voraussetzungen 


Die Blade Element Theory beruht auf einigen Voraussetzungen. Diese Voraussetzungen 


werden im folgenden zunächst aufgelistet und anschließend beurteilt. 


T- . 
Voraussetzungen: 


E: 


1 


I 


Der Druck pı weit vor dem Windrad ist gleich dem Druck p; weit hinter dem 


Windrad. 
Die Strömung ist inkompressibel. 


Der induzierte Widerstand an der Flügelspitze wird nur pauschal für das ganze Rad 
berücksichtigt. 
Es gibt keine Beeinflussung der Strömung an einem Flügelelement infolge von an- 


deren Elementen am gleichen Flügel. 


Es gibt keine Beeinflussung der Strömung an einem Flügelelement infolge von an- 


deren Elementen an anderen Flügeln. 


. Winkeländerungen in der Anströmrichtung durch die Abminderung der Strömungs- 


geschwindigkeit bzw. durch die Aufweitung der Stromröhre sind vernachlässigbar. 


. Die Abminderung der Strömungsgeschwindigkeit beginnt schon vor dem Rotor. Die 


axiale Strömungsgeschwindigkeit in der Rotorebene v, ist der Mittelwert aus den 


Strömungsgeschwindigkeiten weit vor dem Rotor v; und weit hinter dem Rotor v; 


u”) = 


vim 


(vı + v3) 


Die Strömung wird schon vor dem Rotor in Rotation versetzt. Die Umfangsge- 
schwindigkeit relativ zum Erdboden der Strömung im Rotor Au; ist die Hälfte der 


Umfangsgeschwindigkeit Aus weit hinter dem Rotor: 


& 
Auzr = 5 Aus 


Es soll nun untersucht werden, welche Fehler sich aus den Voraussetzungen ergeben: 


zul. Unter der Bezeichnung ”weit vor dem Windrad” versteht man: ”so weit vor dem 


Windrad, daß vom Windrad noch keine Auswirkungen auf die Strömung in diesem 


Punkt meßbar sind”. Man muß etwa mit einer Entfernung des 10fachen Rotorradius 


rechnen. Ein unterschiedlicher Luftdruck kann sich aufgrund von meteorologischen 


Faktoren erst über größere Entfernungen bemerkbar machen und wirkt sich daher 


noch nicht bei Entfernungen von x 10: R aus. Die Voraussetzung kann also als 


gerechtfertigt angesehen werden. Natürlich herrscht kurz vor dem Rotor leichter 


Überdruck und kurz hinter dem Rotor leichter Unterdruck. Da der Kontrollraum 


aber außerhalb dieses engen Bereiches verläuft, braucht die Druckdifferenz nahe am 


Rotor in der Windradberechnung nicht extra behandelt zu werden. Der Druckver- 


lauf in der Umgebung eines Windrades ist in Bild 2.2 gezeigt. 


zu 2. Bis zu einer Machzahl von 0.4 kann mit inkompressibler Strömung gerechnet werden. 


Windräder werden so gebaut, daß eine Machzahl von 0.4 nicht überschritten wird. 


Die höchsten Machzahlen werden an der Flügelspitze erreicht. Bei schnellaufen- 


den Windrädern und hohen Windgeschwindigkeiten kann es ausnahmsweise dazu 


kommen, daß die Forderung Ma < 0.4 nicht mehr erfüllt ist. Näheres kann man [2] 


S.40 entnehmen. 


“ 


zu 3. Der induzierte Widerstand wird durch die Einführung eines Blattzahlverlustfaktors 


berücksichtigt, der sich nicht auf ein Flügelelement bezieht, sondern pauschal für das 


ganze Windrad gilt. Glauert [9] hat jedoch auch eine Rechenmethode angewandt, 


bei der die Wirkung des induzierten Widerstandes auf jedes Flügelelement berück- 


sichtigt wird. Dieses Rechenverfahren wird hier nicht weiter verfolgt. Die Fragen 


zum Blattzahlverlustfaktor werden in einem eigenen Abschnitt noch eingehender un- 


tersucht. Festzuhalten ist hier, daß der induzierte Widerstand nur näherungsweise 


in die Rechnung eingehen kann. 


zu4. Die Strömung an einem Flügelelement wird sicher durch die Strömung an einem 


ku | 


zu d. 


zu 6. 


anderen Flügelelement am gleichen Flügel beeinflußt werden. Z.B..ist es denk- 
bar, daß sich eine noch anliegende Strömung an einer Stelle eines Flügels leichter 
ablöst, wenn die Strömung an einem benachbarten Flügelelement schon abgelöst 
ist. Unterschiedliche Abminderungen der Windgeschwindigkeit werden zu radialen 
Ausgleichsströmungen führen. Vorgänge dieser Art sind in einem einfachen Rechen- 


modell jedoch nicht zu erfassen. 


Bei hohen Drehzahlen eines Windrades mit vielen Flügeln großer Flügeltiefe wer- 
den die Flügel in den Nachlauf des vorauslaufenden Flügels und somit in eine schon 
beeinflußte Strömung gelangen. Auch dieser Effekt ist in einem einfachen Modell 
nicht zu erfassen. Andererseits wird ein Windrad mit niedriger Drehzahl, wenigen 
Flügeln und kleiner Flügeltiefe gar nicht in der Lage sein, die ganze durch die Wind- 
radfläche durchtretende Luft zu beeinflussen. Verluste dieser Art werden durch den 


Blattzahlverlustfaktor näherungsweise berücksichtigt. 


Wie in Bild 2.3 gezeigt, tritt die Luft nicht senkrecht durch die Rotorfläche. Die 


daraus resultierenden Auswirkungen sind gering. 


Die Voraussetzung Nr.7 besagt, daß die Windgeschwindigkeit vor dem Rotor die 
gleiche Abminderung erfährt wie hinter dem Rotor. Mit der einfachen Impuls Theo- 
rie nach Betz kann gezeigt werden, daß die Voraussetzung Nr.7 erfüllt ist. (Siehe 
dazu [8] S.10). Betz gibt jädsch an, daß dieses Ergebnis nur dann gilt, wenn infolge 
der Krümmung der Strombahn nicht zu hohe Zentrifugalkräfte auftreten, d.h. wenn 
die Schnellaufzahl A = NR/vı groß genug ist. Glauert ([9] S.191) hat dies näher 


untersucht. Er definiert eine Gesamtabminderung 


U — Vz 
be 


v1 
und eine Teilabminderung für den Bereich vor dem Rotor 


U — 0 
Wr —— 
v 


zu 8. 


Für den Fall, daß Voraussetzung Nr.7 exakt erfüllt ist, ergäbe sich damit 


ee 
b 


Glauert zeigt jedoch, daß für den Zusammenhang zwischen a und b gilt: 
1-.a)b 
ER. k ei en 


Der zweite Term in der Klammer ist normalerweise klein, so daß es erlaubt ist, hier 


näherungsweise b = 2a zu setzen. Damit erhält man 


1 
als zzetne 


Ersetzen wir nun noch die Teilabminderung a durch den später ausschließlich ge- 


brauchten Abminderungsfaktor € = v3 /vı, so ist 





Tabelle 2.1 zeigt einige Werte dieser Funktion. Man erkennt, daß für Schnellaufzah- 
len A = NR/vı, die größer als 2 sind, Voraussetzung Nr.7 in guter Näherung erfüllt 
ist. 

Angemerkt werden muß noch, daß Glauert dieses Ergebnis erhielt, indem er für 
die Drallströmung einen Potentialwirbel voraussetzte. Das Problem ist, daß das 
Strömungsfeld eigentlich nicht von vornherein festgelegt werden darf, da es erst aus 
einer Superposition einer großen Anzahl von Wirbelfäden entsteht, die sich alle ge- 
genseitig beeinflussen. Es bleibt hier jedoch wohl nichts anderes übrig, als zunächst 


einmal mit dieser Annahme die Kontrollrechnung durchzuführen. 


Voraussetzung Nr.8 kann wieder aus der einfachen Impuls Theorie hergeleitet wer- 
den. Dies hat Hütter gezeigt. ([5] S.18). Da nun der Drall aus der Abminderung 
berechnet werden kann (G1.8), ist die Voraussetzung Nr.8 automatisch erfüllt, sofern 


Voraussetzung Nr.7 gilt. 


2.4 Grundlagen 


Die Formeln zur Blade Element Theory werden für jedes Radelement aus den Grundglei- 
chungen der Strömungsmechanik gewonnen: Kontinuitätsgleichung, Energieerhaltungssatz 
und Impulssatz. Es handelt sich um eine eindimensionale Betrachtung einer inkompres- 
siblen, reibungsfreien Strömung. Mit Hilfe der empirisch ermittelten Daten Gleitzahl e 
und Auftriebsbeiwert c, können die Verhältnisse zweidimensionaler Strömung an einem 
Flügel berücksichtigt werden. Durch die Gleitzahl und den Auftriebsbeiwert geht eine 
Aussage über Zähigkeit, Turbulenzgrad und Reynoldszahl in die Gleichungen ein, die das 
Verhalten des Windrades beschreiben. Die induzierte Wirbelbildung, die vom Druckunter- 
schied an der Profilober- und -unterseite an den Flügelenden herrührt, wird durch einen 
gesonderten Wirkungsgrad, den Blattzahlverlustfaktor n; berücksichtigt. Es werden damit 
die Ergebnisse der Windradberechnung der tatsächlichen dreidimensionalen Strömung an- 
gepaßt. Für jedes Radelement müssen zwei Wirkungsgrade berechnet werden: der ideale 
Leistungsbeiwert cp; und der Flügelwirkungsgrad np. Es wird über allen Radelemente der 
Windradfläche A ein Mittelwert dieser Wirkungsgrade gebildet, so daß man einen neuen 


Wirkungsgrad n erhält: 


n 


neh cepilijnp(i) AA(t) ; (2.2) 


ı=1 


1 


Der Gesamtleistungsbeiwert cp des Windrades hat die Definitionsgleichung 





Antriebsleistung 
e = 
= Gesamtleistung in der Strömung ’ (Def.2.1) 
Lz 
cp=— 
Pi 


Die Gesamtleistung in der Strömung P% ist 


1 
P=> pr) A k (2.2) 


Die Windradleistung steigt also mit der dritten Potenz der Windgeschwindigkeit. A ist die 
gesamte Fläche des Windrades. Bei einem Windrad, das im Nabenbereich keine (profilier- 


ten) Flügel aufweist und dem Wind im Nabenbereich demzufolge keine Energie entziehen 
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kann. muß also schon definitionsgemäß mit einem niedrigeren Wirkungsgrad gerechnet 


werden. Der Gesamtleistungsbeiwert wird in der Blade Element Theory mit 
[PZN (2-3) 
berechnet und man erhält die Windradleistung P aus Def. 2.1 mit 


Pie, (2.4) 


In Bild 2.4 ist die Strömung an einem Flügelelement mit den dort auftretenden Geschwin- 
digkeitskomponenten gezeigt. Bild 2.5 erklärt die Winkelbezeichnungen an einem Flügel- 


element: Windangriffswinkel a, Blatteinstellwinkel 3 und Windeinfallswinkel p. 


Es ist üblich einige dimensionslose Kennzahlen einzuführen. Die in dieser Arbeit verwen- 


deten Bezeichnungen wurden in Anlehnung an Hütter [11] gewählt. 


Abminderungsfaktor E=v/o (Def. 2.2) 
Drallfaktor o=us/u (Def. 2.3) 
Schnellaufzahl IA=u/u (Def. 2.4) 
effektive Schnellaufzahl I, = ef; (Def. 2.5) 
Gleitzahl e= Calw (Def. 2.6) 


Der Index 1 steht immer für die Ebene weit vor dem Windrad. Entsprechend bedeuten 
2 und 3: Windradebene bzw. Ebene weit hinter dem Windrad. Die axiale Geschwindig- 
keit vı ist die Geschwindigkeit der unbeeinflußten Strömung: v; = vo. Bei konstanter 
Anströmgeschwindigkeit kann man in der Schnellaufzahl A eine dimensionslose Drehzahl 


sehen. 
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2.5 Der ideale Leistungsbeiwert cp; 


Der schon erwähnte ideale Leistungsbeiwert cp; ist wie folgt definiert: 





entnommene Leistung aus der Strömung (Def.2.7) 
“L i ef.2.7 


cp; = 
3 gesamte Leistung in der Strömung 


Der ideale Leistungsbeiwert wird nach der Formel 
Berti tee-1) (2.5) 


berechnet. Eine Herleitung dieser Gleichung findet man in [10] S.20. Für einen Horizon- 


talachsrotor ist die Schnellaufzahl A am Flügelelement 


A=\:- (2.6) 


= 
R 
Mit A, ist die Schnellaufzahl an der Flügelspitze (tip) gemeint. Für einen Rotor mit 
geneigter Achse und Flügeln, die gegebenenfalls nicht senkrecht zur Drehachse stehen, ist 


(1:5 — sinrt cos 6) Vcot?s +1 


cosr — sint -cotk :cos® 


A= 





Diese Gleichung ist in [15] S.9 hergeleitet. Die Winkelbezeichnungen Achswinkel r, Ko- 
nuswinkel x und Umlaufwinkel @ sind in Bild 2.6 dargestellt. 


Da die relative Umfangsgeschwindigkeit in der Rotorebene uz unbekannt ist, kann der 


Drallfaktor nicht aus seiner Definitionsgleichung berechnet werden, sondern wird aus 


o=Y1+(1-&)/ (2.8) 


berechnet. Herleitung: [10] S.20. 


Bild 2.7 zeigt, daß der höchste ideale Leistungsbeiwert den Wert 0.59 hat, der bei höheren 
Schnellaufzahlen und einem Abminderungsfaktor € = 1/3 erreicht wird. Für Schnellauf- 
zahlen größer als etwa A = 3 besteht kaum noch eine Abhängigkeit des idealen Leistungs- 
beiwertes von der Schnellaufzahl, sondern fast nur noch eine Abhängigkeit vom Abmin- 


derungsfaktor. Bild 2.8 zeigt den gleichen Sachverhalt noch einmal, jetzt aber mit der 
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Schnellaufzahl als Parameter. Aus Bild 2.7 erkennt man, daß sich bei niedrigen Schnellauf- 
zahlen das Optimum des Abminderungsfaktors vom Wert 1/3 auf den Wert 1/2 verschiebt. 


Zur Berechnung des optimalen Abminderungsfaktors gibt Gourieres ([16] 5.88) die Formel 
1 
RN Et cos( 7 arctanA+ n) (2.9) 


an. 


Warum können nun Windräder mit niedrigen Schnellaufzahlen dem Wind nicht so viel 
Leistung entziehen? Die Leistung P eines Windrades errechnet sich ganz allgemein aus 
dem Drehmoment M: 


[bj 


P=M-u=M-Z X 


Bei gleicher Windgeschwindigkeit vı und gleichem Radius R muß also ein Langsamläufer 
gegenüber einem Schnelläufer ein größeres Drehmoment aufweisen. Nach dem Drehim- 
pulserhaltungssatz übt das Windrad auf die Strömung das gleiche Drehmoment aus, wie 
umgekehrt die Strömung auf das Windrad. Bei einem Langsamläufer verläßt die Strömung 
das Windrad demnach mit starkem Drall, d.h. großer Rotationsenergie, die für eine Nut- 


zung verloren ist. 
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2.6 Der Flügelwirkungsgrad np 


Nicht alle Energie, die dem Wind entzogen wurde, kann zum Antrieb (z.B. eines Genera- 
tors) genutzt werden, weil die Luft den Flügeln bei ihrer Drehbewegung einen Widerstand 
entgegensetzt. Zunächst wird nur der Widerstand am einzelnen Flügelelement berücksich- 


tigt. Die Definitionsgleichung des Flügelwirkungsgrades lautet 


Antriebsleistur.g bei unendlicher Flügelstreckung (Def2.3) 
ef.2. 





ie entnommene Leistung aus der Strömung 


Der Flügelwirkungsgrad wird aus 


2\(e -— Aa) 


Sroa+e) = 


berechnet. Diese Gleichung wurde in [10] S.23 hergeleitet. 


Die Schnellaufzahl in der Rotorebene A, kann nicht direkt aus der Definitionsgleichung 
berechnet werden, weil die Geschwindigkeit u nicht bekannt ist. Man berechnet A, nach 


[10] S.23 aus 
l+o 
1+£ 
Bild 2.9 zeigt den Flügelwirkungsgrad als Funktion der Schnellaufzahl A mit der Profil- 


ı=A- (2.11) 





gleitzahl e als Parameter. Für die Berechnung dieser Kurven wurde für jede Schnellaufzahl 
der optimale N £opı nach Gl. 2.9 gewählt. Aus Bild 2.9 ersieht man, 
daß für Windräder mit hohen Schnellaufzahlen gut profilierte Flügel mit einer hohen Gleit- 
zahl unerläßlich sind. Bei Langsamläufern mit einer Schnellaufzahl von z.B. A; = 2 reicht 
hingegen ein einfaches Profil (eine gewölbte Platte, e = 20) durchaus. Die Leistungsaus- 
beute könnte in diesem Fall auch mit einem aufwendigen Profil nicht wesentlich gesteigert 


werden. 
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2.7 Berechnung des Abminderungsfaktors £ 


Bisher wurde noch nicht behandelt, wie der Abminderungsfaktor € ermittelt wird. Der 
Abminderungsfaktor kann nicht aus seiner Definitionsgleichung berechnet werden, weil die 
Axialgeschwindigkeit hinter dem Rotor vz unbekannt ist. Der Abminderungsfaktor muß 
aus der Rotorgeometrie berechnet werden, denn die Gestalt des Rotors ist es letztlich, die 
die durchströmende Luft in bestimmter Art beeinflußt. Unter Berücksichtigung von Drall 


und der Widerstandskraft am Flügel errechnet sich der Abminderungsfaktor zu 


g=1 - SR. /RHI (+0) x (2.12) 


dnr 
Dabei ist t die Flügeltiefe 
= die Flügelzahl 
r der Radius des Flügelelementes. 


Gl. 2.12 ist im Anhang hergeleitet. Für e — © entspricht diese Gleichung denen aus der 


Literatur ([10] S.33-35,[16] S.88), die für einen Flügel ohne Widerstand hergeleitet wurden. 


1 


Bei Berechnung des Abminderungsfaktors £ muß folgendes beachtet werden: 
e Der Auftriebsbeiwert c„ und die Gleitzahl e sind vom Windangriffswinkel «a abhän- 


gig: 


a = arctan(1/A2) — B ; (2.12a) 


e Die effektive Schnellaufzahl A, ist abhängig von der Schnellaufzahl A und wiederum 


von dem Abminderungsfaktor € (Gl. 2.11). 


Eine Lösung für den Abminderungsfaktor £ ist also nur iterativ möglich. 





U Um die Auswirkungen des Terms l/e in Gl.2.12 zu ermitteln, wurde die Leistung des Windrades 
"Adler 25” der Firma Köster (D-2240 Heide) berechnet. Für die Leistung des Windrades wurde bei 
Berücksichtigung des Term 1/e ein um weniger als 1% geringerer Wert erhalten als beie — 00. Der 


Term kann also vernachlässigt werden, ohne daß ein nennenswerter Fehler auftritt. 
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2.83 Bestimmung des Flügelumrisses 


Aus Gl.12 erkennt man, daß man einen kleinen Abminderungsfaktor (bzw. eine große 
Abminderung der Strömungsgeschwindigkeit) erhält, wenn die Flügeltiefe, die Flügelzahl 
und der Auftriebsbeiwert groß sind. Um an jedem Radius r einen optimalen Abminde- 
rungsfaktor &,pı zu erhalten, muß die Flügeltiefe im Nabenbereich größer sein als an der 
Flügelspitze. Bild 2.10 zeigt den optimalen Flügelumriß im Vergleich zu zwei leichter zu 
fertigenden Kompromißformen. Die Form des Trapezflügels ist eine gute Näherung zum 
optimalen Flügelumriß. Durch die Vorgabe einer bestimmten Tiefenverteilung muß der 
Auftriebsbeiwert in Gl. 2.12 so gewählt werden, daß wiederum der optimale Abminde- 
rungsfaktor erreicht wird. Dies geschieht durch einen veränderten Anströmwinkel a oder 
durch Profile im Nabenbereich, die einen höheren Auftriebsbeiwert haben. Durch diese 
Maßnahmen vergrößert sich natürlich der Profilwiderstand und somit verringert sich die 
Gleitzahl. Damit wird aber der Flügelwirkungsgrad geringer und letztlich die Leistung 
des Windrades. Nach [14] vermindert sich die Leistung bei einem Trapezflügel jedoch nur 
um etwa 1 %. Bei einem Rechteckflügel ist mit Verlusten im Bereich von 7 % zu rechnen. 
Bei der Auslegung eines Trapezflügels soll man so vorgehen, daß die Flügeltiefe an der 
Blattspitze £; der optimalen Flügeltiefe entspricht, während man am Radius r = 0.2R den 


Wert ty2 =1.3...1.7-t, für die dortige Flügeltiefe wählt. [16]. 
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2.9 Bestimmung des Blattwinkels 3 


An der Flügelspitze, wo die optimale Flügeltiefe gewählt werden kann, wählt man auch den 
Auftriebsbeiwert c, durch den Windangriffswinkel a so, daß sich die beste Profilgleitzahl e 
ergibt. Der Profilwirkungsgrad np wird dadurch günstig beeinflußt. Bei dieser Einstellung 
des Profils wird sich ein Auftriebsbeiwert ergeben, der vom maximalen Auftriebsbeiwert um 
einige Anstellwinkelgrade entfernt liegt. Dies ist wichtig, damit bei plötzlicher Zunahme 
der axialen Anströmgeschwindigkeit in einer Böe die Strömung am Profil nicht abreißt. Das 
Abreißen der Strömung bei jeder kleinen Böe würde zu beträchtlichen Leistungseinbußen 
führen, weil die Strömung am Profil eine gewisse Zeit benötigt, um sich wieder an das 
Profil anzulegen. Bei sehr starken Böen ist das Abreißen der Strömung jedoch erwünscht, 
weil dadurch die vom Wind aufgenommene Leistung begrenzt wird und die Anlage vor 


Beschädigungen geschützt werden kann. ' 


Während an der Flügelspitze optimale Bedingungen gewählt werden können, wird bei ei- 
nem vereinfachten Flügelumriß in der Form des Trapezes die Flügeltiefe zur Nabe hin 
immer stärker von der Idealform abweichen. Der aus Gl. 2.12 mit € = £opı gefundene 
Auftriebsbeiwert muß nun durch eine Veränderung des Windangriffswinkels a erreicht wer- 
den. Mit der sich entlang des Rotorblattes in radialer Richtung verändernden Zuströmung 


ändert sich auch der Windeinfallswinkel p. Der Blattwinkel 
B=ypy-a (2.13) 
muß diesen Bedingungen angepaßt werden: 
ß = arctan(1/X2) — a(c.) (2.14) 


Da der Auftriebsbeiwert aus Gl. 2.12 nicht nur vom Abminderungsfaktor und der Geome- 


trie abhängt, sondern auch von der Gleitzahl, muß c. iterativ gefunden werden. 


l Bei dieser vor allem von dänischen Anlagen mit Erfolg eingesetzten Regelung ist zu beachten, daß 
nur die Leistung vermindert werden kann, nicht aber der Schub auf das Windrad. Windräder, die so 


geregelt sind, können zwar auf eine Flügelverstellung verzichten, müssen aber sehr robust gebaut sein. 
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2.10 Der Blattzahlverlustfaktor n: 


Während der ideale Leistungsbeiwert cp; und der Flügelwirkungsgrad np für ein einzelnes 
Radelement gelten, ist der Blattzahlverlustfaktor n. ein Wirkungsgrad, der pauschal die 
Verluste durch den induzierten Widerstand an den Flügelspitzen berücksichtigt. Der Blatt- 
zahlverlustfaktor bezieht sich daher auf das ganze Windrad. Der Blattzahlverlustfaktor 


hat die Definition 


Antriebsleistung 


(Def.2.9) 





aaa Antriebsleistung bei unendlicher Flügelstreckung 


In der Literatur sind verschiedene Formeln zur Berechnung des Blattzahlverlustfaktors 
enthalten: Betz [8] gibt an, daß ac zur Berücksichtigung des rdtebrien Widerstandes 
für das Windrad einen reduzierten Durchmesser annehmen müßte. Rechnet man diese 
Reduzierung des Durchmessers in einen Wirkungsgrad um, so erhält man 


2 

0.44 1 

n: — 1l-= 2 e . (2.15) 
2 y1+%% 


Hütter [11] verwendet einen ähnlichen Blattzahlverlustfaktor: 





2 
1.39 I 
„it l.= —. = s - j 2.16) 
. | z %Jl+ =) \ 
Der Unterschied zur Gl. 2.15 liegt im wesentlichen darin, daß hier die Schnellaufzahl statt 


der effektiven Schnellaufzahl herangezogen wird. Gouriere hat in [16] 


2 

1.394 ä 

„=1i1=- (2.17) 
z  /1+ 2.247-.X2, 


angegeben. Hier wird der Blattzahlverlustfaktor also wieder über die effektive Schnellauf- 








zahl berechnet, die jedoch mit einem Faktor versehen ist. Hennchen berechnet einen Blatt- 


zahlverlustfaktor von 
Lam im Ar 2e 


; 2.1 
2rn-z:-R 2e- X, a8) 





n„=1- 


Der Index m kennzeichnet hier über dem Radius gemittelte Werte. 
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Der Blattzahlverlustfaktor nach Hütter ist in Bild 2.11 dargestellt. Geringe Verluste stellen 
sich demnach bei hohen Schnellaufzahlen und hoher Flügelzahl ein. Bei hoher Flügelzahl 
ergeben sich nach Gl. 2.12 entsprechend schmale Flügel mit großer F lügelstreckung, für 
die der induzierte Widerstand gering ist. Bei unendlich hoher Schnellaufzahl bewegen sich 
bei optimaler Flügelgestaltung unendlich schmale Flügel über die Rotorfläche, die in der 


Lage sind. die Energie des Windes optimal ”abzuernten”. 


Ausnahmsweise ganz anders, als bisher dargelegt, geht Glauert [9| vor. (Wilson und Lissa- 
man [13] haben den Berechnungsgang von Glauert in ihrem Programm zur Windradberech- 
nung angewandt.) Es wird der Einfluß des induzierten Widerstandes auf jedes Radelement 
berücksichtigt. Der induzierte Widerstand hat demnach zur Folge, daß sich Abminde- 
rung und Drall vermindern. Statt einer axialen Differenzgeschwindigkeit v; — v3 erhält 
man nur noch eine axiale Differenzgeschwindigkeit F(vı — v3) und statt einer tangentia- 
len Differenzgeschwindigkeit us — u, erhält man F(us — u,). Dabei ist F der sogenannte 
F-factor oder tip loss factor mit F < 1. Mit anderen Worten erhält man also einen neuen 


Abminderungsfaktor £err und einen neuen Drallfaktor o.yyr mit 
Esel Fi) 


Ger =1+ F(o -1) 


Der tip loss factor F ist 


2 = er 
F = -arccose ? Rein» 
Tr 


Der tip loss factor ist eine Funktion der effektiven Schnellaufzahl (p = arctan(1/A2) ), 
der Flügelzahl = und der radialen Position des Radelementes r. F nimmt einen niedrigen 
Wert an (d.h. hohe Verluste) bei kleiner Flügelzahl, niedriger effektiver Schnellaufzahl und 
einer Position des Radelementes nahe an der Flügelspitze. Bei der Anwendung des tip loss 
factor’ s ist zu beachten, daß der Windeinfallswinkel p von der effektiven Schnellaufzahl 
A, abhängt, die aus A, £ und a berechnet wird. (Gl. 2.11). Der tip loss factor kann also 


nur iterativ zusammen mit anderen Größen der Windradberechnung ermittelt werden. 
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2.11 Der Gesamtlei:tungsbeiwert cp 


Wenn wir nun die drei Wirkungsgrade idealer Leistungsbeiwert cp;, Flügelwirkungsgrad 
7» und Blattzahlverlustfaktor 7. gemeinsam betrachten, so erhalten wir ein teilweise wi- 
dersprüchliches Bild. Hohe Werte für den idealen Leistungsbeiwert und für den Blatt- 
zahlverlustfaktor ergeben sich bei großen Schnellaufzahlen. Demgegenüber gebietet der 
Flügelwirkungsgrad, eher bescheidene Schnellaufzahlen zu wählen. Die Windenergie wird 
heute vor allem zur Elektrizitätserzeugung eingesetzt. Die Generatoren benötigen hohe 
Drehzahlen. Da große Getriebe mit hoher Übersetzung teuer sind und weil im Extrem- 
fall Probleme mit Selbsthemmung im Getriebe auftreten können, ist man an einer hohen 
Drehzahl des Windrades interessiert und somit an einer hohen Schnellaufzahl. Auch das 
Bestreben, größere Windräder zu bauen, führt zu höheren Schnellaufzahlen; denn da 
w-R 


v 


A= 





ist, steigt bei gleicher Drehzahl die Schnellaufzahl proportional zum Radius des Windra- 
des. Aus diesen Gründen ist heute eindeutig eine Tendenz zu größeren Schnellaufzahlen 
festzustellen. Um einen guten Flügelwirkungsgrad zu erreichen, muß mit den besten Profi- 
len gearbeitet werden, die hohe Gleitzahlen aufweisen. Diese aufwändig profilierten Flügel 
sind teuer. Die Kosten können jedoch durch eine geringere Flügelzahl reduziert werden. 
Dieser Zielsetzung kommt auch entgegen, daß Flügel ohnehin aus Festigkeitsgründen nicht 
zu schmal und dünn sein dürfen und daher nach Gl. 2.12 zur Erzielung eines optimalen 


Abminderungsfaktors &.p: eine geringere Flügelzahl nötig ist. 


Eine theoretische Grenze der Entwicklung größerer Windräder mit höheren Schnellaufzah- 
len ist durch die Machzahl mit dem Wert 0.4- Ma gegeben. So liegt danach z.B. für einen 
Rotor bei einer Windgeschwindigkeit von 7 m/s! die höchste Schnellaufzahl bei A; = 19.6. 
Beim GROWIAN wurde als höchste Anströmgeschwindigkeit der Wert 0.29- Ma festgelegt. 
Bei einer Auslegungswindgeschwindigkeit von 11 m/s ergab sich danach eine Schnellaufzahl 


14=9. [14,19]. 








l Mittleres Jahresmittel der Windgeschwindigkeit auf Sylt nach [18]. 
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2.12 Der Momentenbeiwert cy 


Der Momentenbeiwert cu hat die Definitionsgleichung 
M 
m Def.2.10 
MTTge?.R-A .. 


und kann aus dem Gesamtleistungsbeiwert über 





cp 
= — 2.19 
MT (2.19) 
berechnet werden. Der zum Anlaufmoment M, gehörende Momentenbeiwert cm, ist 
BER: 3% 
we IN, für dı=0 f (2.20) 


Windräder mit hoher Schnellaufzahl haben ein kleines Drehmoment. Schwierigkeiten ver- 
ursacht bei ihnen das ebenfalls kleine Anlaufmoment, welches durch eine für den Anlauf 
falsche Einstellung der Flügel verursacht wird. Schnellaufende Windräder müssen daher 
entweder durch Fremdenergie angeworfen werden oder eine Verstellmöglichkeit der Flügel 
aufweisen. Die Strömung am Flügel bei verschiedenen Betriebszuständen zeigt Bild 2.12. 
Ein weiterer Grund für das schlechte Anlaufverhalten von Schnelläufern ist die geringe 


Rotorblattflächendichte 





(2.21) 


wobei Ar die Summe der in die Rotorebene projezierten Blattflächen ist. 


2.13 Der Widerstandsbeiwert cw 


Der Widerstandsbeiwert cyr eines Windrades ist nicht zu verwechseln mit dem Wider- 


standsbeiwert c„, eines Flügels. Der Widerstandsbeiwert cw eines Windrades wird durch 


Foa 


cw = —— (Def.2.11) 
zpuf A 
definiert und kann aus dem Abminderungsfaktor berechnet werden: 
1x R 
ew=7 2 (1-£?(r,6)) AA . (2.22) 


(Herleitung in|3]). Die Axialkraft F,. ist eine wichtige Größe für die Auslegung des Turmes. 
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2.14 Überschlagsrechnungen 


Bei Überschlagsrechnungen kann man vom Verhalten eines Radelementes auf das ganze 
Rad schließen. Dazu teilt man sinnvollerweise den vom Flügel bestrichenen Kreisring, 
der vom Radius R, bis zur Flügelspitze mit dem Radius R geht, in zwei flächengleiche 


Kreisringe. Diese Kreisringe berühren sich beim Radius 
AT ı 2 2 22 
lg (R? + RZ) ; (2.23) 


welcher als effektiver Radius für das eine Radelemement der Überschlagsrechnung gewählt 
wird. 

Hütter [10] gibt einen effekiven Radius r.yr = 0.72R ohne Beweis an. Zu bemängeln ist an 
dieser Formel, daß sie nicht den Radius der Flügelwurzel R, berücksichtigt. In der Formel 
von Hütter ist ein Wert R„/R = 0.19 vorgegeben, welcher wohl den meisten Windrädern 


entspricht. 


Bei Rechnungen mit optimalen Werten, basierend auf dem optimalen Abminderungsfaktor 


£opt, kann Tabelle 2.2 hilfreich sein. 


DD 
[SS 


3 Beschreibung des Panelverfahrens 


3.1 Allgemeine Beschreibung 


Das Panelverfahren hat seinen Namen vom englischen Wort ”panel”, was in diesem Zu- 
sammenhang so viel heißt wie Plättchen oder Verkleidungsblech. Das Panelverfahren ist 
ein umfangreiches Programmsystem und dient zur Berechnung der Strömung um beliebige 
Körper. Das Programm ermöglicht ebenfalls, die Strömung um mehrere (am Institut für 
Mechanik (IM) maximal zwei) Körper mit den auftretenden Interferenzerscheinungen zu 
berechnen. Die Oberfläche der zu berechnenden Körper wird dabei näherungsweise als aus 
ebenen Plättchen (panels) zusammengesetzt angenommen. Bild 3.1 zeigt dies am Beispiel 
einer Flügel-Rumpf-Kombination. Aus Bild 3.1 erkennt man, daß die Panels an den Stellen 
kleiner gewählt werden, an denen die Geometrie des Körpers oder der Verlauf der zu er- 
wartenden Strömung sich stark ändern. Je größer die Anzahl der Panels ist, desto genauer 
kann die tatsächliche Form des Körpers wiedergegeben werden, womit auch die Genauig- 
keit der Lösung zunimmt. Es ist jedoch zu beachten, daß die Rechenzeit quadratisch mit 


der Anzahl der Panels steigt. [20]. 


Das Panelverfahren ermöglicht die Berechnung der räumlichen Potentialströmung. Ähn- 
lich wie in der einfachen Potentialtheorie, wo umströmte Körper durch eine Überlage- 
rung der Elemente: Parallelströmung, Quelle, Senke, Potentialwirbel und Dipol modelliert 
werden, geschieht dies auch beim Panelverfahren. Für die Oberflächenpanels wird eine 
konstante Quell-Senkenverteilung angenommen, deren Quellstärke so bestimmt wird, daß 
senkrecht zum Panel kein Materiestrom fließt. Die Überlagerung (Superposition) der kon- 
stanten Quell-Senkenverteilung aller Panels zusammen mit der Anströmung ergibt dann 


die Strömung um den Körper. 


Erst durch den Einsatz leistungsfähiger elektronischer Datenverarbeitungsanlagen konn- 
ten die ersten Panelverfahren in den 60er Jahren realisiert werden. Zunächst war es nur 
möglich, Verdrängungsprobleme zu lösen (z.B. Flugzeugrumpf), später wurden dann auch 


Auftriebsprobleme (Flügel) behandelt. 


Um einen Flügel mit Auftrieb zu berechnen, muß neben der Oberfläche auch noch die 
Abflußströmung vorgegeben werden. Aus der Theorie heraus kann die Strömung an der 
Profilhinterkante nicht bestimmt werden. Bild 3.2 zeigt dazu drei verschiedene Möglichkei- 
ten: In Bild 3.2a wird die Hinterkante von der Unterseite her umströmt, in Bild 3.2b von 
der Oberseite her. Bei derartigen Umströmungen einer scharfen Kante entstehen große 
Strömungsgeschwindigkeiten und Unterdrücke. Wie die Erfahrung zeigt, tritt am Flügel 
in erster Näherung die Strömung aus Bild 3.2c auf. Die Strömung fließt an der Hinterkante 


glatt ab. ! Diese Tatsache wird Kuttasche Abflußbedingung genannt. 


Um die Kuttasche Abflußbedingung zu berücksichtigen, werden sogenannte Wirbelleiter 
modelliert, die sich aus numerischen Gründen in das Tuneze des Kirosie hinein fortsetzen. 
Bild 3.3 zeigt die zum Körper aus Bild 3.1 gehörende "innere Fläche”. Die innere Fläche 
ist bei einem Flügel die Skelettfläche und wird ähnlich der Oberfläche in ”innere” Panels 
eingeteilt. Außerhalb des Flügels werden die Wirbelleiter bis etwa zu einer Entfernung der 
5-fachen Flügeltiefe fortgesetzt. Die Steigung der äußeren Wirbelleiter ist die Steigung der 
Tangente an der Skelettfläche in der Profilhinterkante. Die Panels der Wirbelleiter sind mit 
Dipolen belegt. Die Dipolstärke eines Panels ist konstant. Die Kuttasche Abflußbedingung 
wird in Kontrollpunkten der Wirbelleiter erfüllt, die aus numerischen Gründen etwa 1% 


der örtlichen Flügeltiefe hinter der Hinterkante liegen. 


Das am Institut für Mechanik vorhandene Programm ”IfM-Panel” hat einen Aufbau, der 


sich an einer von MBB entwickelten Programmversion orientiert.[20]. 


! Durch die Zähigkeit wandert der Staupunkt etwas nach oben. Der Auftrieb wird dadurch etwas 


vermindert. 


3.2 Hauptgleichungen 


Für den Benutzer des Panelverfahrens ist das Programm "IfM-Panel” im wesentlichen eine 
"black box”. Eine vollständige Beschreibung des Programms ist am Institut vorhanden. 
Eine Übersicht über die Hintergründe des Panelverfahrens kann man [21] und [22] entneh- 
men. Die Hintergründe des Panelverfahrens sind hier im Hinblick auf die Besonderheiten 


der Windradanwendung wiedergegeben. (Siehe dazu auch [23]). 


Der Windradflügel ist ein Körper, der sich um eine Achse dreht. Für einen Beobachter, der 
vor dem Windrad auf der Erde steht, ist dieser Vorgang nicht stationär. Grund: Der Flügel 
bewegt sich am Beobachter vorbei und bietet demnach ein mit der Zeit sich änderndes Bild. 
Der Vorgang läßt sich dadurch ukänr machen, daß ein Koordinatensystem gewählt wird, 
welches mit dem Flügel fest verbunden ist. Vorausgesetzt werden soll zunächst, daß die 


Anströmung des Windrades im Absolutsystem üp homogen und achsparallel ist: 
Up = (Vo ‚0, 0) 


Beim Übergang vom Absolutsystem ins Relativsystem wird der Parallelanströmung vp eine 
Geschwindigkeit überlagert, die der Geschwindigkeitsverteilung eines rotierenden Festkör- 
pers entspricht. Diese Strömungsform wird Solid-Body-Strömung genannt. Das Strö- 


mungsfeld der Solid-Body-Strömung ist 


UR=-w-ixr 


Dabei ist w die Winkelgeschwindigkeit des Windrades 
i der Einheitsvektor in x-Richtung 
r der Vektor zu irgendeinem Punkt auf dem Windrad. 


Das Vektorfeld der (unbeeinflußten) Anströmung im Relativsystem vw ergibt sich aus der 


vektoriellen Addition von Parallelströmung vp und Solid-Body-Strömung vr: 


vw =tp+Vr 
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Das Vektorfeld der tatsächlichen Strömung um den Windradflügel im Relativsystem V 
setzt sich zusammen aus dem Beitrag der Anströmung vUw und dem eines überlagerten 


Störfeldes grady, welches durch den umströmten Körper verursacht wird: 


_ 


V=vw + gradp : (3.1) 


Die kinematische Strömungsbedingung [21] auf der Oberfläche des Rotors verlangt, daß 


die Strömung nicht in die Oberfläche eindringt, d.h. 
V-a=0 . (3.2) 


Dabei ist n der Normalenvektor auf einem Flächenelement dA der Flügeloberfläche. Mul- 


tipliziert man Gl. 3.1 mit n und berücksichtigt Gl. 3.2, so ist 
0=vw-n + grad(p)-n 


oder 


oder 


grad(p)-n=-vp-n-vür:n ; (3:3) 


Ausgangspunkt der dreidimensionalen, inkompressiblen, stationären Strömungsberechnung 


mit der Potentialtheorie ist die Laplace-Gleichung 


RE 








a 7 Bar ee 
die aus der Definition des Potentials 
Op dp dp 
u-= J v= — » w= — 
dr dy 2 
und aus der Kontinuitätsgleichung 
du x dv dw 0 
e: u 


26 


folgt. Zu beachten ist, daß die Laplace-Gleichung nur auf das Störpotential % angewen- 
det werden kann. Das Potential des Vektorfeldes V existiert nicht, weil es aufgrund der 


überlagerten Solid-Body-Strömung nicht wirbelfrei ist. x 


Das Störpotential einer räumlichen Quell-Senkenverteilung für einen beliebigen Punkt P 


pas(P) = —/° ers, (3-4) 


Dabei ist Sx die Körperoberfläche und dSx ein Flächenelement davon, 


beträgt nach [21 


r der Abstand von P zu einem beliebigen Punkt auf $x und 
co die Quell-Senkenstärke eines Punktes auf Sk. 


Dieses Störpotential genügt zur Darstellung eines reinen Verdrängungseffektes. Um einen 
Körper mit Auftrieb zu berechnen, sind ferner die Dipole der Wirbelleiter zu berücksich- 
tigen, deren Stärke aus der Abflußbedingung zu bestimmen ist. Das Störpotential eines 


räumlichen Dipols für einen beliebigen Punkt P ist nach [21] 


ep(P) 2a (Sn), 2 ) dSn (3.5) 


Dabei ist Sn die Oberfläche des Nachlaufes und dS n ein Flächenelement davon 
r der Abstand von P zu einem beliebigen Punkt auf Sx 
u die Dipolstärke eines Punktes auf Sn 
8/ön der Differentialquotient in Richtung der Normalen von Sn. 


Gl. 3.4 und Gl. 3.5 sind Lösungen der Laplace-Gleichung. Diese Lösungen erhält man 


nach Bronstein ([26], 5.553) durch Anwendung der Greenschen Funktion. 





I Soll die Strömung im Relativsystem direkt berechnen werden, so geht dies nur über die Poisson- 
Gleichung ([24], S.103) 

ar 00 U 

O2? Oy? 8:2 


2w . 





wobei P die Stromfunktion ist. 


[S°) 
| 


Das Gesamtstörpotential ergibt sich aus der Überlagerung von Gl. 3.4 und Gl. 3.5 zu 


Eingesetzt in (3.3) errechnet sich für Punkte auf der Körperoberfläche Q 


grad |- R Ile un ax + [| m (Sn) (2) asx| so) - 


= -tp(Q)-R(Q) - er(Q)-R(Q) (3.6) 





Dies ist eine Deliflineare Integralgleichung für die Singularitätenstärken o($x) und u(Sn). 
Da die Oberfläche des Körpers aus einzelnen Panels mit konstanter Singularitätenbelegung 


besteht, lassen sich die Doppelintegrale vereinfachen: 


/). as = De, 8 45; 
Wurst en [l,&0)“ 


Dabei bezeichnet m die Anzahl der Körperpanels und n die Anzahl der Panels für den 


Nachlauf. Führt man nun noch einige Abkürzungen ein 


p 

ee 

Fo Ale 
Kj 

EN = NO 

ET AV 


Agj = [-uraa /f. "as, .n(Q) 
re), (5) | 


-tp(Q) n(Q) - Er(Q)  R(Q) 
Vo 








Ro = 


und teilt (3.6) durch V., um die Geschwindigkeit zu normieren, so erhält man, durch 


Einführen der Abkürzungen in (3.6) nach dem Umordnen die vereinfachte Schreibweise 


YX,-Ay+) YBo=Re - (3.7) 


3-1 3-1 
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Dies ist ein lineares algebraisches Gleichungssystem mit (m + n) unbekannten Dipol- bzw. 
Quell-Senkenstärken. Der kritische Punkt des Programms zum Panelverfahren ist die 
Auflösung dieses umfangreichen Gleichungssystems. Beim Gauß-Algorithmus steigt der 
Rechenaufwand mit der dritten Potenz der Anzahl der Unbekannten. Dieses Verfahren 
ist daher nur schlecht anwendbar. Ausgenutzt werden kann, daß die Diagonalelemente der 
Matrix von (3.7) dominant sind. Dies ermöglicht die Anwendung der Gauß-Seidel-Methode, 
bei der der Lösungsvektor iterativ berechnet wird [27]. Je mehr die Diagonalelemente in 
ihrer Größenordnung überwiegen, desto schneller konvergiert das Verfahren. Nach [21] läßt 


sich das Konvergenzverhalten positiv beeinflussen, wenn 
- die Winkeländerung von der Ebene eines Panels zur Ebene des Nachbarpanels klein 
ist 
- die Änderung der Panelabmessungen von einem Element zum anderen klein ist. 


Eine entscheidende Reduzierung der Rechenzeit ergibt sich vor allem durch das Ausnutzen 
von Symmetriebedingungen. Beim Windrad ist dies die Rotationssymmetrie. Es braucht 


dann nur die Strömung um einen Flügel berechnet zu werden. 


Nachdem das Gleichungssystem gelöst ist, ist das gesamte Störpotential bekannt. Es 
können die Geschwindigkeiten an der Körperoberfäche Vx berechnet werden und daraus 


nach Bernoulli der Differenzdruck 
(3.8) 


Der Auftriebsbeiwert ergibt sich dann aus einer Summation der einzelnen, für jedes Panel 


berechneten Kraftkomponente. 


3.3 Aufbau des Programms 


Mit dem Programm ’IfM-Panel” kann die Strömung um maximal zwei auftriebsbehaftete 
Körper berechnet werden. Die Geometrie der Körper muß vom Benutzer in Unterprogram- 
men berechnet oder eingelesen werden. Die Geometrie des ersten umströmten Körpers wird 
über das Unterprogramm KOERPI eingegeben, die Geometrie des zweiten Körpers (falls 
vorhanden) über das Unterprogramm KOERP2. Die Anzahl der Koordinaten, die für die 
Panels eingegeben werden können, ist auf 2500 begrenzt. Die Felder, die die Koordinaten 
enthalten sind mit (50,50) dimensioniert. Die x-, y- bzw. z-Koordinaten eines Paneleck- 
punktes werden in den Feldern PX(I,J), PY(I,J) bzw. PZ(I,J) gespeichert. Die Lage des 
Eckpunktes im Gitter wird durch: die I- und J-Koordinaten beschrieben, wie aus Bild 3.4 


ersichtlich ist. Ein Panel wird also von vier benachbarten Eckpunkten gebildet. 


Die Geometrie der Wirbelleiter muß bei Körpern mit Auftrieb über die Unterprogramme 
WLEIT1 und WLEIT?2 bereitgestellt werden. WLEIT1 enthält die Wirbelleiter für den 
ersten Körper und WLEIT2 gegebenenfalls die Wirbelleiter für den zweiten Körper. Die 
Eckpunkte der Wirbelleiterpanels werden in den Feldern PXW(1,J), PYW(1,J), PZW(1L,J) 
gespeichert. Der Index I kennzeichnet die Nummer der Sprosse, der Index J kennzeichnet 
die Nummer des Holmes. Die Holme weisen in Strömungsrichtung, also in Richtung der 
Profiltiefe, während die Sprossen dazu senkrecht stehen und in Richtung der Flügellänge 
weisen. 

Steuerparameter, Daten der Parallelanströmung und Geometriedaten werden vom Pro- 
gramm von einer Datei gelesen. In diese Datei muß der Benutzer auch die Daten für sein 
spezielles Problem eingeben, die von den Unterprogrammen KOERP und WLEIT ange- 
fordert werden. Es ist möglich, der parallelen Anströmung, andere Elemente wie Quellen 
und Dipole zu überlagern oder das Relativsystem rotieren zu lassen. Der Einfluß dieser 
Singularitäten auf die Geschwindigkeitskomponenten in den Panel-Nullpunkten wird in 


einem Unterprogramm ADDSIN berechnet. 


Das Programm berechnet die Druckverteilung auf der Körperoberfläche, den Auftriebs- 
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beiwert und den Widerstandsbeiwert. Die aus dem umfangreichen Gleichungssystem be- 
rechneten Singularitätenstärken werden für spätere Auswertungsläufe abgespeichert. Dies 
geschieht ebenfalls mit Daten, die zur graphischen Weiterverarbeitung erforderlich sind. 


Mit Programmen, die zur Datenaufbereitung dienen, sind verschiedene Ausgaben möglich: 
- Plotten der Druckverteilung 
- Plotten von Stromlinien 
- Plotten des Fernfeldes der Strömung 


- Darstellung der Strömung durch Vektorpfeile. 
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4 Unterprogramme zum Panelverfahren 


4.1 Das Unterprogramm KOERP1 zur Berechnung der Flügelpanels 


Aufgabe des Unterprogramms KOERPI ist es, aus Geometriedaten, die von einer Da- 
tei eingelesen werden, die Paneleckpunkt-Koordinaten zu berechnen. Die x-, y- und z- 


Koordinaten der Paneleckpunkte werden in den Feldern PX(1,J), PY(I,J) bzw. PZ(1,J) 


abgespeichert. 


Als erstes muß ein Koordinatensystem festgelegt werden. Das hier gewählte - mit dem 
Windrad fest verbundene - Koordinatensystem zeigt Bild 4.1. Die Parallelanströmung 
erfolgt längs der x-Achse. Der Windradflügel dreht sich beim Blick in Windrichtung gegen 


den Uhrzeigersinn. 


Ein Windradflügel besteht in den meisten Fällen nicht nur aus einem Flügelprofil. Um den 
unterschiedlichen Anströmbedingungen zwischen der Flügelwurzel und der Flügelspitze ge- 
recht zu werden, finden unterschiedliche Profile Verwendung. Reine Profile aus dem Pro- 
filkatalog liegen jeweils nur in einem Querschnitt vor. Querschnitte des W indradflügels zwi- 
schen zwei reinen Profilen sind Übergangsprofie, die die reinen Profile möglichst gleichmäßig 
ineinander überleiten sollen. Das Unterprogramm KOERP1 liest die benötigte Information 
über die Profile aus einer Eingabedatei ein. Die Daten werden in diese Datei so eingegeben, 
wie sie aus Profilhandbüchern entnommen werden: Alle Profildaten sind auf die Profiltiefe 
t bezogen. Die Skelettlinie des Profils entspricht der x-Achse mit dem Wert z = O in der 
Profilnase und dem Wert z = 1 in der Profilhinterkante. Die Profilsaugseite wird durch 
Wertepaare (20,yo) beschrieben und die Profildruckseite durch Wertepaare (zu,yu). Die 
Lage der reinen Profilquerschnitte auf dem Windradflügel wird durch den Radius r zur 
Windradachse angegeben. 

Die Flügeltiefe t wird durch Wertepaare (r,t) beschrieben. In der gleichen Art erfolgt dies 
für die Flügelverwindung mit den Wertepaaren (r,A3). (Siehe auch Bild 2.5). Aß ist 
der Verwindungswinkel in bezug zum Winkel an der Flügelspitze (r = R), der den Wert 
Aß(R) = 0 hat. Die Winkelstellung des gesamten Flügels wird durch den Winkel der 
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Blattspitze 3, angegeben. Es ist also: 
B(r) = ı + AA(r) . (4.1) 


Bei Windrädern mit Flügelverstellung kann die Winkelstellung des gesamten Flügels wäh- 
rend des Betriebs variiert und dadurch an die variablen Größen Windgeschwindigkeit, 


Drehzahl und Lastmoment angepaßt werden. 


Sch'ießlich werden noch der Windradradius R und der Radius der Flügelwurzel R„ be- 
nötigt. Die Werte r, R, Ru, und t sind dimensionsbehaftet. Sie können in beliebiger 
Weise eingegeben werden, müssen jedoch maßstäblich zueinander sein. Die Werte werden 
im Unterprogramm KOERP1 dimensionslos gemacht, indem sie auf den Windradradius R 


bezogen werden. 


Da es nicht unbedingt sinnvoll ist, sämtliche Panels gleich groß zu wählen (siehe dazu 
Kapitel 3.1 und 3.2), kann der Benutzer eine Panelgrundverteilung angeben. Nach dem 
Studium einiger Panelverteilungen aus [28] wird als eine mögliche Verteilung der Panels in 


Richtung der Flügeltiefe eine Grundverteilung 
0.06, 0.14, 0.30, 0.30, 0.20 


vorgeschlagen, die in Bild 4.2 gezeigt ist. Für eine Verteilung der Panels in Flügellängs- 


richtung erscheint folgende Grundverteilung sinnvoll: 
0.06, 0.09, 0.11, 0.12, 0.12, 0.12, 0.12, 0.11, 0.09, 0.06 


Eine derartige Flügelaufteilung ist in Bild 4.1 am Beispiel des Windrades Nr. 7 der Ae- 
rodynamischen Versuchsanstalt Göttingen gezeigt. Die hier vorgeschlagenen Panelgrund- 
verteilungen sollen dem Anwender lediglich als Hilfestellung dienen und können je nach 
Problemstellung anders gewählt werden. Das Unterprogramm VERTEIL berechnet zu ei- 
ner gewählten Grundverteilung und der gewünschten Panelzahl eine Panelendverteilung. 
Die oben angegebene Grundverteilung in Richtung der Flügeltiefe besteht z.B. aus 5 Wer- 
ten; bei einer gewünschten Panelzahl in dieser Richtung von 15 wird jedes Element der 


Grundverteilung noch einmal in 3 Panels unterteilt. 
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Nachdem jetzt die geforderte Verteilung der Panels auf dem Flügel festliegt, geht es darum, 
die Paneleckpunkt-Koordinaten zu berechnen. Zunächst müssen aus den Profildaten die 
Werte yo und yır berechnet werden, die zu den x-Werten der Panelendverteilung gehören. 
Dies wird mittels Spline-Interpolation durchgeführt. Ein Unterprogramm SMOSPL exi- 
stiert dazu am IfM, welches in einem Umdruck [29] beschrieben ist. Das Unterprogramm 
ermöglicht eine glättende Spline-Interpolation. Bei der hier auftretenden Problemstel- 
lung geht es jedoch nicht um eine Glättung der Werte, sondern darum, daß die Spline- 
Funktionen direkt durch die gegebenen Punkte der Profilkontur verlaufen sollen. Wird der 
Glättungsparameter in SMOSPL auf den Wert $ = 1.0!" gesetzt, so führt das Unterpro- 


gramm keine Glättung der Funktionswerte durch. 


Da die Paneleckpunkte in der Regel nicht dort liegen, wo sich die reinen Profilquerschnitte 
des Flügels befinden, müssen Werte der Übergangsprofile berechnet werden. Auch dies 
geschieht mit Hilfe des Unterprogramms SMOSPL. In der Praxis kommt es häufig vor, 
daß die Profile direkt an der Flügelspitze und an der Flügelwurzel nicht bekannt sind. Da 
mit der Spline-Interpolation nur Zwischenwerte berechnet werden können, werden Über- 
gangsprofile in den Flügelaußenbereichen durch lineare Extrapolation im Unterprogramm 


ESPLINE berechnet. 


Bevor die Paneleckpunkt-Koordinaten schließlich den Feldern PX, PY, PZ zugewiesen 
werden können, muß noch eine Koordinatentransformation durchgeführt werden, in der 
der Verlauf der Flügeltiefe und die Flügelverwindung berücksichtigt werden. Die Koor- 
dinatentransformation übernimmt das Unterprogramm KOORD. Die in diesem Unter- 
programm benötigten Werte t(r) und AA(r) werden ebenfalls durch Extrapolation und 


Spline-Interpolation aus den Eingabedaten gewonnen. 
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4.2 Das Unterprogramm WLEIT1 zur Berechnung der Wirbelleiter 


Da der Windradflügel ein Körper mit Auftrieb ist, muß ein Nachlaufmodell erarbeitet wer- 
den. Zur Erfüllung der Kutta-Joukowski Bedingung an der Profilhinterkante muß die Lage 
der Wirbelleiter und der Kontrollpunkte berechnet werden. Die x-, y- und z-Koordinaten 
der Paneleckpunkte der Wirbelleiterpanels werden in den Feldern PXW(LJ), PYW(1LJ) 
bzw. PZW(I.J) abgespeichert. 


Ein ideales Nachlaufmodell (Wirbelleitermodell) muß so aussehen, daß die Wirbelleiter 
genau den Stromlinien folgt, die im Staupunkt an der Profilhinterkante beginnen. Es 
tritt hier das Problem auf, daß das Strömungsfeld dazu als bekannt vorausgesetzt werden 
müßte, obwohl es doch gerade erst darum geht, das Strömungsfeld zu berechnen. Für ge- 
naue Rechnungen kann nur iterativ vorgegangen werden: Für den ersten Iterationsschritt 
wird ein einfaches Nachlaufmodell angenommen. Für die Berechnung des nächsten (ge- 
naueren) Nachlaufmodells dienen die Werte des Strömungsfeldes aus der vorhergehenden 
Rechnung. In dieser Arbeit wird keine Iteration des Nachlaufmodells durchgeführt. Das 


Programm WLEIT1 ermöglicht die Wahl zwischen zwei Wirbelleitermodellen: 


1. Wirbelleitermodell: 


Zur Berechnung der Wirbelleiter wird angenommen, daß Luftteilchen, die nach der Um- 
strömung des Flügels die Profilhinterkante verlassen, ihre Umfangs- und Axialgeschwin- 
digkeit beibehalten. Der Quotient aus Umfangs- und Axialgeschwindigkeit dieser Teil- 
chen wird der Steigung der Skelettlinie in der Profilhinterkante gleichgesetzt. Es wird 
angenommen, daß die Umfangsgeschwindigkeit proportional zum Radius ist und daß die 
Radialgeschwindigkeit Null ist. Die angenommenen Stromlinien, denen die Wirbelleiter 
folgt, haben demnach einen konstanten Abstand zur Drehachse des Windrades. Vom Pro- 


grammanwender wird der Umfangswinkel Ymaz festgelegt, über den sich die Wirbelleiter 
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erstrecken soll. sowie die Anzahl der Panels der Wirbelleiter. 


Die genannten Annahmen zur Konstruktion der Wirbelleiter sind in der Praxis nicht erfüllt. 
Bild 2.4 zeigt: Die Axialgeschwindigkeit nimmt ab, die Umfangsgeschwindigkeit nimmt 


zu und der Radius der Stromlinien (im Bild ist eine Stromlinie gestrichelt gezeichnet) 


vergrößert sich. 


Bei einer Flügelverwindung 


Aß(r) = arctan(1/Az(r)) 


erhält man als Wirbelleiter eine Helix. die man am sinusförmigen Verlauf erkennt. Bild 4.3 
und Bild 4.4 zeigen dazu je ein Beispiel. In Bild 4.4 hat der Flügel die Flügeltiefe t = 0. 
Derartig geordnete Verhältnisse ergeben sich in den meisten Fällen nicht. Bild 4.5 zeigt 
perspektivisch ein Beispiel einer Wirbelleiter, bei der der Bereich der Wirbelleiter an der 


Flügelspitze dem Nabenbereich vorauseilt. 


2. Wirbelleitermodell: 


Das 2. Wirbelleitermodell basiert auf dem ersten. Auch bei diesem Modell wird das 
Verhältnis aus Umfangs- und Axialgeschwindigkeit für den Bereich direkt hinter der Profil- 
hinterkante (Q,) au: der Steigung der Skelettlinie in der Profilhinterkante berechnet. Wäh- 
rend jedoch beim 1. Wirbelleitermodell angenommen wird, daß das Verhältnis Umfangs- 
zu Axialgeschwindigkeit (Q) längs des Verlaufs der Wirbelleiter den Wert der Skelettlinien- 
steigung Qı beibehält, wird beim 2. Wirbelleitermodell angenommen, daß das Verhältnis 
Umfangs- zu Axialgeschwindigkeit den Wert annimmt, den ein optimal gestaltetes Wind- 
rad in der Rotorebene aufweist: Q2. Das Verhältnis von Umfangs- zu Axialgeschwindigkeit 
in der Rotorebene ist nach Def. 2.5 u2/v2 = A2 = Q:- A, kann für ein optimal 
gestaltetes Windrad aus Gl. 2.11 mit Gl. 2.9 und Gl. 2.8 berechnet werden. Für den 


Übergang des Wertes Q} auf den Wert Qz längs des Drehwinkels der Wirbelleiter y wird 
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ein linearer Verlauf angenommen. Um zu gewährleisten, daß die Wirbelleiter nicht in ei- 
nen anderen Windradflügel hineinläuft, wurde festgesetzt, daß nach einem Drehwinkel der 
Wirbelleiter >. = 360°/(2-=) der Wert Qı in den Wert Qz übergegangen sein muß. 
Es ist demnach: 


Q: ‚= Qı er 
er für TE% 


Te 
Qaly) = Q: für 12Y%e 
Auch dieses Wirbelleitermodell entspricht noch nicht den in Bild 2.4 gezeigten Verhältnis- 


sen. Die Nachlaufmodelle können nur durch (richtige) Ergebnisse gerechtfertigt werden. 


4.3 Berechnung von Kennzahlen 


Im Programm ”IfM-Panel” wird ein Druckbeiwert c, berechnet, der auf die achsparallele 


Anströmgeschwindigkeit V„ bezogen ist: 
el- i (4.2) 


Die Geschwindigkeit tangential zur Körperoberfläche Vx ist, verursacht durch die Rota- 
tion des Flügels, an der Flügelspitze deutlich größer als Vo. Demzufolge ergeben sich dort 
betragsmäßig hohe Druckbeiwerte. In der Propellerberechnung bezieht man den Druck- 
beiwert daher auf die Geschwindigkeit V2, + r? -w? (vergleiche mit [30]) und erhält einen 
modifizierten Druckbeiwert 


V2 
ed K 


PO v2 + r.02 


In einem Unterprogramm CPNEU wird der Druckbeiwert c, in den modifizierten Druck- 
beiwert c; nach der Formel 


v2 


2 2.w2 
vn + rw 


f=1+(%,-]1) (4.4) 
umgerechnet. Durch die Einführung des modifizierten Druckbeiwertes ist es möglich, 
Druckbeiwerte des Windradflügels aus der 3dimensionalen Rechnung mit 2dimensionalen 


Profilberechnungen zu vergleichen. 


Im Unterprogramm CAFLE wird der Auftriebsbeiwert eines Flügelelementes berechnet. 
Unter einem Flügelelement soll hier ein Flügelabschnitt verstanden werden, der aus einer 
Reihe von Panels besteht, die alle die gleich z-Koordinate (Bild 4.1) haben. Windradachse 
haben. Die Auftriebsbeiwerte der Flügelelemente werden aus den modifizierten Druckbei- 
werten berechnet. Der Auftriebsbeiwert bezieht sich auf die unbeeinflußte Anströmung des 


Flügels vw. 
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Im normierten Windrad, so wie es im Programm vorliegt, gelten folgende Werte: 


RE Bi ; 
R=>18 ; 
A=r » 
A, =w 


Diese Werte wurden den Berechnungen der Windradkenngrößen zugrunde gelegt. Die Bei- 
werte für die Flügelelemente beziehen sich auf die Fläche der zugehörigen Radelemente 


(vergleiche mit Bild 2.1). Beiwerte für das ganze Windrad beziehen sich auf die Wind- 
radfläche A= r. 
Der Widerstandsbeiwert cw (des Windrades, nicht des Flügelprofils) wird aus der Defini- 


tionsgleichung (Def. 2.11) berechnet: 


Ip) Art) at). (4.5) 


cw = 


A| 


Dabei ist z die Flügelzahl 

i eine Laufvariable 
m die Anzahl der Körperpanels 
cp(i) der Druckbeiwert am i-ten Panel 
Ap(i) die Panelfläche des i-ten Panels 
ri(i) der Normalenvektor auf der Panelfläche des i-ten Panels 
i der Einheitsvektor in x-Richtung. 

Der Widerstandsbeiwert des Flügelelementes cw,s wird ähnlich berechnet: 


zZ 
- 


cwE = Ar(r) 2 


ı 
>" -%(i)- Apt)-nl) ei . (4.6) 


ı=1 





Dabei ist Ag(r) die Fläche des Radelementes mit dem mittleren Radius r 
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! die Anzahl der Panels eines Flügelelementes. 


Aus den Widerstandsbeiwerten wird jeweils ein Abminderungsfaktor £ berechnet. Der 


Abminderungsfaktor stellt eine integrale Größe über die Fläche dar, für die der Wider- 


standsbeiwert berechnet wurde. 


Der Momentenbeiwert cm wird aus der Definitionsgleichung (Def. 2.10) berechnet 


ah, 


CM 


u c>(i)- Apli)- Ri) -3 


Dabei ist j der Einheitsvektor in y-Richtung. 


Der Momentenbeiwert für das Flügelelement cm,g ist 


I 
OELOEzEE (4.8) 





CM,E = 
i=1 


Aus dem Momentenbeiwert wird nach Gl. 2.19 der Leistungsbeiwert cp des Windrades 


und des Flügelelementes cp,& berechnet. 
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5 Ergebnisse des Panelverfahrens 


5.1 2dimensiorale Profilberechnung mit dem Panelverfahren 


Zunächst ging es darum, zu prüfen, welche Genauigkeit bei der Anwendung des Panelver- 
fahrens zu erwarten ist. Zu diesem Zweck wurde die Strömung um Flügelprofile berechnet 
und mit Literaturangaben verglichen. In der Literatur sind Auftriebs-, Widerstandsbei- 
wert und Druckverteilung bei 2dimensionaler Umströmung von Flügelprofilen angegeben. 
Eine 2dimensionale Flügelumströmung wurde in der Rechnung mit dem Panelverfahren 


dadurch modelliert, daß eine Flügellänge von der 100fachen Flügeltiefe gewählt wurde. 


Bild 5.1 zeigt die berechnete Druckverteilung um ein Profil NACA 2412 bei drei verschiede- 
nen Anstellwinkeln «. In Bild 5.2 sind die entsprechenden gemessenen Druckverteilungen 
gezeigt. Für Anstellwinkel von —1.7° und 7.4° besteht eine sehr gute qualitative und 
quantitative Übereinstimmung von Messung und Rechnung; lediglich im Bereich der Pro- 


filhinterkante bestehen leichte Abweichungen im Druckverlauf. 


Wie Bild 5.2 zeigt, ist die Strömung bei einem Anstellwinkel von 17.9° über einem weiten 
Bereich der Flügeltiefe abgelöst. Mit dem Panelverfahren wird jedoch eine anliegende 
Strömung berechnet, wodurch sich auf der Profilsaugseite erheblich höhere Unterdrücke 
ergeben. Der Druckverlauf an der Profilhinterkante stimmt nicht mit dem gemessenen 


überein. 


Die Strömung an dem symmetrischen Profil NACA 0015 wurde für die Anstellwinkel 0°, 5° 
und 10° berechnet. Die Druckverteilungen zu diesen Anstellwinkeln zeigt Bild 5.3. Über 
die Auftriebsbeiwerte bei diesen Anstellwinkeln ist ein Vergleich der berechneten Werte 
mit den gemessenen Werten aus der Literatur möglich (Bild 5.4). Da bei 10° die Strömung 
offensichtlich noch nicht abgelöst ist, besteht für alle drei Anstellwinkel eine sehr gute 
Übereinstimmung der berechneten Werte mit den Meßwerten. Für 0° und 5° konnte selbst 
die Strömung an der Profilhinterkante recht genau vorausberechnet werden. Bei 10° weist 


der Druckverlauf in der Profilhinterkante Ungenauigkeiten auf. 
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Bei allen Rechnungen ergaben sich für die beiden Panels auf der Profildruck- und -saugseite 
an der Profilhinterkante die gleichen Druckbeiwerte. Obwohl der berechnete Druckbeiwert 
an dieser Stelle nicht immer dem gemessenen entsprach, so zeigt der jeweils gleiche Wert 
doch an, daß der Abströmvorgang durch die gewählte Wirbelleiter (siehe 3.1) gut modelliert 
wurde. 

Der Widerstandsbeiwert eines Profils wird mit dem Panelverfahren prinzipiell immer zu 
Null berechnet. Abweichungen von dem Wert Null sind äußerst klein und haben numerische 


Gründe. 


Die folgende Liste zeigt eine Zusammenstellung berechneter und gemessener Auftriebsbei- 


werte: 


Profil Ce,» berschnet C„, gemessen 
NACA 2412 0.024 
NACA 2412 0.862 
NACA 2412 0.950 
NACA 0015 0.0 
NACA 0015 0,50 
NACA 0015 0.35 





Die gemessenen Auftriebsbeiwerte für das Profil NACA 2412 sind [31] entnommen. Sie 
gelten für eine Reynoldszahl von 2.7 - 10%. Die Auftriebsbeiwerte für das Profil NACA 
0015 entstammen [12]. Sie wurden bei einer Turbulenz von 2% und einer Reynoldszahl 
von 6.6-10° gemessen. Für kleine Anstellwinkel stimmen die Ergebnisse der Rechnung gut 
mit den gemessenen Werten überein. Die Abweichungen sind kleiner als die, die man durch 
unterschiedliche Reynoldszahlen in der Messung erhält. Sobald sich jedoch die Strömung 
ablöst, ist der Auftriebsbeiwert aus der Rechnung viel zu hoch und daher nicht mehr zu 


verwerten. 


42 


5.2 Beschreibung der untersuchten Windräder 


Im Rahmen dieser Arbeit wurden zwei Windräder mit dem Panelverfahren nachgerechnet: 


Das Windrad Adler 25 


Das Windrad Adler 25 wird von der Maschinenfabrik Köster gebaut. Entwickelt wurde 
es von der DFVLR unter dem Namen DEBRA 25 mit Mitteln des BMFT. Der Rotor des 
Windrades ist in Lee des Turmes angeordnet. Das Windrad hat drei Flügel mit einer 
Schränkung von 8°, einen Rotordurchmesser von 25 m und eine Nennleistung von 100 kW. 
133]. Aus Konstruktionsunterlagen [34] der Firma Köster geht hervor, daß sechs verschie- 
dene Flügelprofile eingesetzt wurden: NACA 4415, NACA 4418, NACA 4421, FX 77 W 
258, FX 77 W 343 und eine modifizierte Form des Profils FX 77 W 405. Die Lage dieser 
Profile auf dem Flügel zeigt Bild 5.5. Für die Rechnung mit dem Panelverfahren wurde der 
Flügel erst ab dem Radius R, = 3.22 m modelliert. Das modifizierte Profil FX 77 W 405 
und der sich zur Nabe hin anschließende Kreisquerschnitt blieben unberücksichtigt. Der 
verbleibende Flügel wurde in 15 Panelreihen unterteilt, die in Bild 5.6 gezeigt sind. Jeweils 
der mittlere Querschnitt der 2., 5., 7., 10. und 15. Panelreihe (von der Flügelwurzel aus 
gezählt) entspricht einem der berücksichtigten Profile. Das Strak des wie in Bild 5.6 pa- 
nelierten Flügels ist in Bild 5.7 gezeigt. Die Profilkoordinaten der NACA-Profile konnten 
mit einem Programm berechnet werden, das am IfM vorhanden ist. Die Profilkoordinaten 


der Profile FX 77 W 258 und FX 77 W 343 sind in [35] veröffentlicht worden. 
Das Windrad Nr. 7 der Aerodynamischen Versuchsanstalt zu Göttingen 
In der IV. Lieferung der Aerodynamischen Versuchsanstalt zu Göttingen [36], die 1932 


erschien, wurden die Ergebnisse der Messung an vier verschiedenen Windrädern veröffent- 
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licht. Das Windrad Nr.7 hat die größte Ähnlichkeit mit modernen Windrädern und wurde 
daher für die Berechnungen im Rahmen dieser Arbeit ausgewählt. Zusammen mit den 
Meßergebnissen ist auch eine Übersichtszeichnung des Windrades veröffentlicht worden 
(Bild 5.8). Das Windrad wurde für Windkanaluntersuchungen gebaut. Es hat daher nur 
einen Durchmesser von 0.9 m. Die zwei Flügel des Windrades weisen eine Schränkung 
von 7° auf. Für die Flügel wurden drei verschiedene Profile eingesetzt: Gö 449, Gö 450 
und Gö 451. Profilkoordinaten für diese Profile konnten nicht ausfindig gemacht werden. 
Eine recht genaue Abbildung der Profile ist jedoch in [37] gegeben. Aus dieser Abbildung 
wurden die Profilkoordinaten abgemessen und anschließend im Programm SMOSPL [29] 
geglättet. Der panelierte Flügel des Windrades wurde bereits in Bild 4.1 gezeigt. Das 


Strak des wie in Bild 4.1 panelierten Flügels zeigt Bild 5.9. 
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5.3 Ergebnisse der Windradberechnug mit dem Panelverfahren 


Ein großer Teil der Programmläufe diente dazu herauszufinden, in welcher Weise der Flügel 
paneliert und die Wirbelleiter modelliert werden muß. Zunächst ging es darum, eine sinn- 
volle Anzahl der Panels in Richtung der Flügeltiefe zu bestimmen. Ausgegangen wurde 
von der Grundverteilung, die in Kap. 4.1 angegeben ist. Bild 5.10 a) zeigt die Druckver- 
teilung am Flügelelement bei nur je 5 Panels auf der Profildruck- und -saugseite, also bei 
nur 10 Panels pro Panelreihe. Die Bilder 5.10 b) und 5.10 c) geben die Druckverteilung 
bei 20 bzw. 40 Panels wieder. Es kann anhand dieser Druckverläufe festgestellt werden, 
daß erst bei einer höheren Panelzahl von 30 bis 40 Panels pro Flügelelement sinnvolle 
Werte berechnet werden. Die Überschneidung der Kurven des Druckverlaufs im Bereich 
der Profilhinterkante verschwindet bei hoher Panelzahl fast vollständig. Der Wert maxi- 
malen Unterdrucks (im Bild nicht mehr gezeigt) beträgt: a) c», = —2.54, b) c, = -3.35, 
c)c» = —4.46. Es ist mit diesen Werten erklärbar, daß die berechnete Windradleistung 


mit steigender Panelzahl pro Flügelelement zunimmt. 


Auch in radialer Richtung können verschiedene Panelzahlen gewählt werden. Bild 5.11 
zeigt die Druckverteilung an einem ausgewählten Radius. Die Druckverteilung wurde im 
Fall a) mit 10 Panels und im Fall b) mit 20 Panels in radialer Richtung berechnet. Für 
beide Fälle ist der ermittelte Druckverlauf identisch. Eine größere Genauigkeit ergibt sich 
im Fall b) bei der Berechnug von Kennzahlen des gesamten Windrades. Um die Kennzahlen 
Widerstandsbeiwert cw, Momentenbeiwert cm und Leistungsbeiwert cp zu erhalten, muß 
nach Gl. 4.5 und Gl. 4.7 der Druckbeiwert c, summiert werden. Am realen Windrad 
ändert sich der Druckbeiwert über dem Radius stetig. In der Rechnung ist dieser stetige 
Verlauf durch eine Treppenfunktion angenähert, da von einem konstanten Druckbeiwert 
auf jedem Panel ausgegangen wird. Durch die Wahl einer hohen Panelzahl in radialer 
Richtung kann die Treppenfunktion an die reale, stetige Funktion angenähert werden, was 


zu einer größeren Genauigkeit bei der Berechnung der Kennzahlen des gesamten Rotors 


führt. 
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Bei der Variation der Geometrie der Wirbelleiter wurden folgende Ergebnisse erzielt: 


Nr. des maximaler Anzahl Gesantleistungss- 
irbelleiter-|Drehwinkel der | der beiwert 
modells Jirbelleiter Sprossen Cp 





Alle Programmläufe wurden für den Kennfeldpunkt X: = 9.0, 3: = —1.2° durchgeführt, für 
den ein Leistungsbeiwert von cp = 0.44 gemessen wurde. Die Ergebnisse verdeutlichen, 
in welch starkem Maße die berechnete Leistung von der Modellierung der Wirbelleiter 
abhängt. Die nach dem 1. Wirbelleitermodell berechnete Wirbelleiter bewegte sich mit 
steigendem Drehwinkel aufgrund des negativen Blattwinkels an der Flügelspitze 3; ge- 
gen (!) die Windrichtung. Das Ergebnis war ein extrem hoher (nach Betz nicht erreichba- 
rer) Leistungsbeiwert. Durch Anwendung des 2. Wirbelleitermodells konnte die Leistung 
nur geringfügig reduziert werden. Ein größerer Einfluß als von der Wahl des Modells geht 
vom maximalen Drehwinkel der Wirbelleiter aus. Ein Drehwinkel von 90° lieferte die 
höchste Leistung. Bei größeren und kleineren Drehwinkeln als 90° wurde eine geringere 
Leistung berechnet. Bei einem Drehwinkel von 720° und 24 Wirbelleitersprossen wurde 
ein Leistungsbeiwert cp = 0.386 berechnet, der niedriger als der gemessene Wert ist. Der 
Drehwinkel von 360° lieferte ein fast genaues Ergebnis. Daß dies nicht zu verallgemeinern 


ist, zeigt Bild 5.12. 


Am Windrad Adler 25 wurde der gemessene Leistungsbeiwert für ausgewählte Punkte, die 
im Diagramm mit O gekennzeichnet sind, nachgerechnet. (Bild 5.12). Alle Programmläufe 
wurden mit dem 2. Wirbelleitermodell durchgeführt. Ergebnisse, die mit einem maximalen 


Drehwinkel der Wirbelleiter Xmaz = 90° und 6 Wirbelleitersprossen berechnet wurden, 
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sind mit ”Y” gekennzeichnet. Ergebnisse der Programmläufe, die mit einem Drehwinkel 
“mar = 360° und 24 Wirbelleitersprossen erzielt wurden. erhielten das Symbol "A” zur 
Kennzeichnung. Die Ergebnisse zeigen. daß eine kurze Wirbelleiter mit Imaz = 90° auf 


keinen Fall ausreichend ist. Auch die Ergebnisse mit der 360°-Wirbelleiter sind noch zu 


ungenau. 


Die Bilder 5.13 bis 5.17 zeigen den Vergleich von 3dimensionaler und ?2dimensionaler Rech- 
nung mit dem Panelverfahren an den 5 äußeren Profilen des Adler 25. Für die 2dimen- 
sionale Vergleichsrechnung wurde aus den Daten der 3dimensionalen Rechnung der An- 
stromwinkel der Profile in der Rotorebene a = aa(r) nach Gl. 2.12a mit Gl. 2.11, Gl. 2.8 
und dem Abminderungsfaktor € aus Gl. 2.22 mit cw = cw,g aus Gl. 4.6 berechnet. Die 
Kurven der 3dimensionalen Druckverteilung weisen alle eine leichte Überschneidung im 
Bereich der Profilhinterkante auf, was darauf hindeuten kann, daß die Wirbelleiter noch 
nicht optimal gewählt wurde. Eine derartige Überschneidung der Kurve des Druckbei- 
wertes tritt in der 2dimensionalen Vergleichsrechnung nicht auf. Die Druckverläufe der 2- 
und 3dimensionalen Rechnung stimmen qualitativ und quantitativ überein, lediglich bei 
den Wortmann-Profilen ergab die 3dimensionale Rechnung eine stärkere Ausprägung des 
Gebietes mit Unterdruck auf der Profildruckseite. (Bilder 5.13 und 5.14). Zur Flügelspitze 
hin nimmt der Unterdruck am Flügel ab. Das hat nach Bild 5.18 auch eine Abnahme des 
Auftriebsbeiwertes zur Folge. Der größere Auftriebsbeiwert in der Nähe der Flügelwurzel 
ist durch einen größeren Anströmwinkel verursacht. Aufgrund der Ausgleichsströmung 
an Flügelspitze und Flügelwurzel und des induzierten Widerstandes wird hier ein Auf- 
triebsbeiwert c„ = 0 erwartet. Bild 5.18 zeigt ein deutliches Abfallen des Auftriebsbei- 
wertes an den Flügelenden, jedoch wird wegen der nicht unendlich kleinen Panelgröße der 
Wert c„. = 0 nicht erreicht. Die Auftriebsverteilung ändert sich umgekehrt proportional 
zum Radius. Eine derartige Auftriebsverteilung entspricht Gl. 2.12 mit Gl. 2.6 und 


€ = Eopt = 1/3 = konst. nach Tab. 2.2. 
Mit dem Panelverfahren wurden zwei Programmläufe mit veränderter Geometrie des Wind- 
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rades Adler 25 durchgeführt. Es wurde der Kennfeldpunkt A; = 6.2, 3; = 6° gewählt. Die 
Wirbelleiter. die nach dem 2. Wirbelleitermodell berechnet wurde. hatte einen Drehwinkel 
von {max = 90° und 6 Wirbelleitersprossen. Bei einer Reduzierung der Flügeltiefe auf die 
Hälfte des ursprünglichen Wertes fiel der berechnete Leistungsbeiwert von cp = 0.23 auf 
cp = 0.13. Eine Leistungsverminderung entspricht der Erwartung, denn nach GI. 2.12 
wird die durch das Windrad strömende Luft bei kleiner Flügeltiefe nicht mehr so stark 


verzögert. 


Im zweiten Programmlauf wurde die Flügelzahl von 3 auf 1 reduziert. Um jedoch in 
erster Näherung die gleiche Verzögerung der Strömung zu erreichen, wurde gleichzeitig 
die 3fache Flügeltiefe gewählt. Der Leistungsbeiwert fiel dennoch vom Wert cp = 0.23 
auf den Wert cp = 0.15. Auch dieses Ergebnis entspricht den Erwartungen, weil durch 
die geringere Flügelstreckung und dem dadurch zunehmenden induzierten Widerstand ein 


Leistungsverlust auftritt. Weitere Parameterstudien sind in Kap. 6 beschrieben. 


Die Ergebnisse der Berechnung des Windrades Nr.7 der Aerodynamischen Versuchsanstalt 
Göttingen zeigt Bild 5.19. Bei der Schnellaufzahl A = 6.0 ist der berechnete Leistungs- 
beiwert cp kleiner als der gemessene. Bei X; = 8.0 wurde ein zu großer Leistungsbeiwert 
berechnet. Der Fehler der beiden berechneten Leistungsbeiwerte ist im Vergleich zu den 
Ergebnissen am Windrad Adler 25 gering. Der berechnete Wert des Widerstandsbeiwertes 


cw ist deutlich niedriger als der gemessene Wert. 


6 Ergebnisse der Windradberechnung nach Hütter 


Ein vorhandenes Programm ”CP” [15], das auf der Windradberechnung nach Hütter [10] 
und [11] basiert, wurde im Rahmen der Diplomarbeit erweitert, so daß aufwendig profilierte 
Windräder wie das Windrad Adler 25 berechnet werden können. Das Programm CP 
ermöglicht auch die Berechnung von Windrädern, die eine nicht horizontale Achse haben. 
Aus konstruktiven Gründen weist das Windrad Adler 25 einen Achswinkel von 10° auf. 
Um eine stationäre Strömung zu erhalten, wurde die Schrägstellung der Achse bei den 
Berechnungen mit dem Panelverfahren vernachlässigt. Programmläufe mit CP zeigten, 
daß der maximale Leistungsbeiwert bei Berücksichtigung von Achs- und Konuswinkel nur 
um weniger als 1% geringer ist als bei einer Vergleichsberechnung des Windrades als reinem 
Horizontalachsrotor. Eine Verminderung des Leistungsbeiwertes macht sich offensichtlich 
erst bei größeren Achswinkeln bemerkbar [15]. Es ist demnach anzunehmen, daß durch 


die Berechnung des Adler 25 als reinem Horizontalachsrotor keine nennenswerten Fehler 


gemacht werden. 


Bild 6.1 zeigt das Ergebnis einer Windradberechnung, wie sie in Kap. 2.4 bis 2.11 be- 


schrieben ist. 


Das Kennfeld des Leistungsbeiwertes cp für das Windrad Adler 25 aus Bild 5.12 wurde 
mit dem Programm CP nachgerechnet. Das Ergebnis der Rechnung zeigt Bild 6.2. Die 
berechnete Leistung ist immer etwas kleiner als die gemessene - nie größer. Für Blattwinkel 
Bı > 10° sind die Abweichungen gering und von der Größenordnung der Genauigkeit der 
Profilkennzahlen Auftriebsbeiwert c, und Gleitzahl e. Für 9, = 14° liegt die berechnete 
Leistung deutlich unter der gemessenen. Auch die berechnete Leistung für kleine Schnel- 
laufzahlen wurde deutlich zu niedrig berechnet. Insgesamt kann jedoch eine recht gute 


Übereinstimmung zwischen Rechnung und Messung festgestellt werden. 


Bild 6.3 zeigt den Einfluß der Flügeltiefe auf den Gesamtleistungsbeiwert cg des Windra- 
des Adler 25. Sowohl bei Verdoppelung als auch bei Halbierung der Flügeltiefe sinkt der 


maximale Leistungsbeiwert. Eine Vergrößerung der Flügeltiefe führt zu einer Verschie- 
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bung des Maximums des Leistungsbeiwertes zu niedrigeren Schnellaufzahlen. Bei einer 


Verkleinerung der Flügeltiefe tritt der umgekehrte Effekt auf. 


Wird das Produkt aus Flügelzahl und Flügeltiefe konstant gehalten (z-t = konst.), dann 
sinkt der Gesamtleistungsbeiwert mit abnehmender Flügelzahl. Die Leistungsverminde- 
rung ist ausschließlich auf den Blattzahlverlustfaktor (Bild 2.11) zurückzuführen. Für den 
idealen Leistungsbeiwert cp; und den Flügelwirkungsgrad np erhält man in jedem dieser 


Fälle den gleichen Wert. 
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7 Berechnung instationärer Strömungsfelder 


Bei einer homogenen und achsparallelen Anströmung eines Windrades erhält man durch 
den Übergang vom Absolutsystem ins Relativsystem im Relativsystem wieder eine sta- 
tionäre Strömung. Die Hauptgleichungen des Panelverfahrens wurden in Kap. 3.2 für 


diesen Fall angegeben. 


Bei Rotoren mit geneigter Drehachse ist jedoch selbst die Strömung im Relativsystem 
noch instationär. Eine instationäre Strömung im Relativsystem erhält man auch, wenn die 
Grenzschicht der Luftströmung über der Erdoberfläche berücksichtigt werlen muß oder bei 
Berechnung der Rotor-Turm- Wechselwirkung. Auch eine Fehlanströmung des Windrades 
aufgrund einer falschen Windrichtungsnachführung der Windkraftanlage verursacht eine 


instationäre Strömung im Relativsystem. 


In diesen Fällen empfiehlt es sich, eine quasi-stationäre Berechnung der Strömung mit dem 
Panelverfahren für verschiedene Umlaufwinkel 8 der Flügel durchzuführen. Die Strömungs- 
geschwindigkeit an der Körperoberfläche Vx wird dann, wie im stationären Fall, nach Gl. 
3.1 berechnet. Ein Unterschied zur stationären Rechnung ergibt sich jedoch bei der Be- 
rechnung des Differenzdruckes (Gl. 3.8) bzw. ee Druckbeiwertes (Gl. 4.2): Während im 
stationären Fall die einfache Bernoullische Gleichung herangezogen werden konnte, muß 
für die instationäre Rechnung die verallgemeinerte Bernoullische Gleichung angewandt 
werden. 
Die Bernoullische Gleichung läßt sich aus der Eulerschen Bewegungsgleichung durch Inte- 
gration längs einer Stromlinie herleiten. Eine einfache Herleitung ist in Stromlinienkoordi- 
naten möglich. Die Integration kann jedoch auch - mit größerem Aufwand - in kartesischen 
Koordinaten erfolgen. Dieser Weg soll hier beschrieben werden. Die Herleitung erfolgt da- 
bei im Hinblick auf die Windradanwendung. Es gilt hier, wie in Kap. 3.2, 
V=öow+Vp (3.1) 


mit 


Uw=tUp+VUR 


ol 


und 
v; R=iw: ix r 
Mit dem Fehlanströmwinkel 9 aufgrund falscher Windrichtungsnachführung, dem Grenz- 


schichteffekt und dem Achswinkel r nach Bild 2.6 gilt für die instationäre Rechnung jedoch 
vp = (Vo: cosd -cosr, Vo - sind, Vo - sınr) 


MtV=(uwi+v.j+w- k) lautet die Eulersche Bewegungsgleichung: 


av ge VL V 
eu u. 





N 
+ FR + gVz=0 . (6.1) 


Mit Vektorschreibweise erhalten wir 





— = (V:Vv)V 


Für die Eulersche Bewegungsgleichung ergibt sich damit 


en -  ÖV 1 
(V-V)V + Er + ut +gVz=0 3 (6.2) 


Während der Zeit dt bewegt sich ein Teilchen in der Strömung auf einer Stromlinie um ds’ 
weiter. Wenn wir nun die Terme aus Gl. 6.2 mit ds’ skalar multiplizieren, so erhalten wir 
eine skalare Gleichung, die die Terme Druck p, Geschwindigkeit V und die Höhe : entlang 


der Stronlinie miteinander verbindet: 
PIERRE BEN e > 
Ba Te De ; (6.3) 
p 
Dabei ist 
d-5=deri + dyj + dzk (längs s) 
Es sollen die Terme aus Gl. 6.3 zunächst einzeln behandelt werden: 


et [3 et 2 [dei + dyj + dzk] 
p 


dx Oy Ö: 
_1[9p Op Op s 
- EZ - y - ni (längs s) 


und 


eu [dei +. dd + dk) 
— gVz-ds= - gdz (längs s) 


Unter Zuhilfenahme einer Rechenregel des Nabla-Kalküls [38] können wir schreiben: 


a $ SE ar e 
(V-VV.d= ver - Vx(Vx 2) . ds 


- 


Der letzte Term auf der rechten Seite dieser Gleichung ist Null, weil V parallel zu ds ist. 


Folglich ist 














an ” RE 
(V-V)V-ds=SV(V-V)-ds (längs) 
1 [-IV?  -ov2  ov2 a 
5 & = + | [tel a + ai] 
1 [3V? öVv? öv? 
| + Bar a 
1 
2 


(V-V)V-.ds=-d(V?) (längs s) 


1 
— ds + — +gd: + zauv?)=0 (längs s) 
pP 
Die Integration dieser Gleichung liefert die verallgemeinerte Bernoullische Gleichung: 


r h 2 . 
ds N 92 + — = konst (längs s) 1 (6.4) 


- 


Im Falle der Windradberechnung beginnt die Integration auf einem Punkt auf der Strom- 
linie im Unendlichen (©) und endet auf einem Punkt der Stromlinie auf der Körper- 


oberfläche (Ä). Betrachten wir zunächst einmal den instationären Term. Mit Gl. 3.1 


ist 





= 
Vo - . ö ” — _ 
= ? Ne sind, sinr) — a + —(Vp)| - ds 
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Eine Vereinfachung dieses Ausdrucks erhält man, wenn sich die Windgeschwindigkeit 
V>(:) nicht mit der Zeit ändert (keine Böen) und wenn das Windrad mit konstanter 


Drehzahl läuft. Es ist dann 


K = K 

OV (0) 

in rn 
ii a /. ee 


K £ _ Ps: PN —_ 
er (+ Ba ER) (de fan i+ del) 


u GM \ De Y z 
He; 
| 


Da die Orte der betrachteten Integrationsgrenzen sich im Relativsystem nicht bewegen, 


= ce ine d 
[ z%-5/ dp = zuloK = Po) 


oo oo 


gilt 


Für das Störpotential ist 90 = 0. Also 


KöV _, Spk 
:dS = —— 


ö dt 


en 
Führen wir nun auch noch die Integration der anderen Terme aus Gl. 6.4 aus, wobei die 
Luftdichte p in der Windradberechnug als konstant angenommen werden kann, so lautet 
das Ergebnis 


ÖpK 1 E 1 e° 
—— + S(V£ - VE) + (Pr -Po)=0 . (6.5) 


Für den Druckbeiwert erhalten wir mit px =p 








_ Vz 
5-8 Peg & 


9} 
Rn: 6.6 
1 p v2 v2 v2 ot ( ) 





Für den Term ?ox muß eine numerische Differentiation durchgeführt werden. Zunächst 
sind in einem bestimmten Intervall des Umlaufwinkels 9 [0 ; 6] Programmläufe 
durchzuführen, wobei der Umlaufwinkel von 6, bis 9 kontinuierlich mit einer Schrittweite 
A® zu erhöhen ist. Bei konstanter Drehzahl entspricht jedem Schritt A® eine Zeitdifferenz 
At. Je nach der Anzahl der Programmläufe können unterschiedlich genaue Formeln zur 
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Berechnung des Terms "5° herangezogen werden: 
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Formel mit einem Fehler 2. Ordnung: 


OPKi _ PKRi+ı — PKi-ı 
öt 2 At 





Formel mit einem Fehler 4. Ordnung: 


dori _-PRis+2 + 8 PrRitı - BpRi-ı + PRi-2 
dt 12 At 


In den Formeln bezeichnet px,; den i-ten Wert der Größe px im betrachteten Intervall. 
In der Literatur sind weitere Formeln zur numerischen Differentiation angegeben, die noch 
kleinere Fehler aufweisen, sowie Formeln für die Berechnung der Ableitung von Funkti- 
onswerten am Rande des betrachteten Intervalls. Die Formeln werden aus den Satz von 


Taylor hergeleitet. Hinweise dazu findet man in [27]. 


8 Ausblick 


Das in dieser Arbeit beschriebene Panelverfahren zur Windradberechnung ermöglicht die 
Ermittlung von Windradkenngrößen, ohne daß dazu Auftriebs- und Widerstandsbeiwerte 
der verwendeten Profile bekannt sein müssen. Die Windradleistung konnte nicht mit der 
erhofften Genauigkeit berechnet werden. Hierzu bedarf es noch einer genaueren Modellie- 
rung der Wirbelleiter oder einer Iteration der Gestalt der Wirbelleiter, wie in Kap. 4.2 
angedeutet. In weiteren Berechnungen könnte versucht werden, die Erkenntnisse über 
den Verlauf der Strömung hinter dem Windrad (wie in Bild 2.4 gezeigt) in die Gestal- 
tung der Wirbelleiter einfließen zu lassen. Weitere Parameterstudien an der Geometrie 
der Wirbelleiter müssen zeigen, ob die Windradleistung gegen einen bestimmten Wert 
konvergiert, wenn die Wirbelleiter immer länger gewählt wird. Die Grenzen der Anwen- 
dung des Panelverfahrens können schon jetzt klar gesehen werden: Genaue Ergebnisse aus 
dem Panelverfahren können nur im Bereich der Auslegungsschnellaufzahl erhalten wer- 
den. Bei zu niedriger Schnellaufzahl wird der Anströmwinkel am Flügel zu groß, es kommt 
dann in der Realität zur Ablösung der Strömung. Da das Panelverfahren immer mit an- 
liegender Strömung rechnet, wird die berechnete Leistung deutlich über der gemessenen 
Leistung liegen. Wie Bild 6.1 zeigt, wird die Leistung bei hohen Schnellaufzahlen durch 
den Flügelwirkungsgrad begrenzt. Das Panelverfahren berücksichtigt nur den induzierten 
Widerstand, nicht aber den Reibungswiderstand, was einem Flügelwirkungsgrad np = 1 
entspricht. Es wird daher erwartet, daß auch bei hohen Schnellaufzahlen die berechnete 
Leistung deutlich über der gemessenen Leistung liegt. Problematisch für die praktische 
Anwendung sind die langen Rechenzeiten des Panelverfahrens. Es wird kaum möglich sein, 
ein ganzes Kennfeld in vertretbarer Zeit zu berechnen. Ein Ausweg könnte hier sein, die 
aerodynamischen Beiwerte der Profile in einer 2dimensionalen Rechnung mit dem Panel- 
verfahren zu ermitteln, um diese dann in einem herkömmlichen Programm zur Windradbe- 
rechnung zu verwerten. Wenn es gelingt. mit dem Panelverfahren unter Berücksichtigung 
der genannten Einschränkungen verwertbare Ergebnisse zu erzielen, ist der Weg offen für 


weitergehende Berechnungen, wie z.B. der Wechselwirkung des Windrades mit anderen 
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Bauteilen. Es können hier Erkenntnisse erlangt werden, die den herkömmlichen, einfachen 


Rechenprogrammen nicht zugänglich sind. 


-1 


or 
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Tabellen 


Tabelle 2.1 


Berechnung des Abminderungsverhältnisses % als Funktion der Schnellaufzahl A und des 


Abminderungsfaktors £ 


A |% 0.5 1.0 2.0 3.0 4.0 5.0 

£=1/3 0.111 0.778 0.944 0.975 0.986 0.991 
{=1/2 0.250 0.813 0.953 0.979 0.988 0.993 
£=2/3 0.444 0.861 0.965 0.985 0.991 0.994 


£=1.0 1.000 1.000 1.000 1.000 1.000 1.000 


[97T] 3uimupsısgpespury op sy sfeuydg 77 >IPqeL 


0.100 0.670 0.473 
0.200 0.768 0.451 
0.300 0.873 0.432 
0.400 0.984 0.416 
0.500 1.099 0.403 
0.600 1.219 0.393 
0.700 1.343 0.384 
0.800 1.470 0.377 
0,900 1.600 0.371 
1.000 1.732 0.366 
1.100 1.866 0.362 
1.200 2.002 0.359 
1.300 2,140 0.356 
1.400 2.279 0.353 
1.500 2.419 0.351 
1.600 2.560 0.349 
1.700 2.702 0.348 
1.800 2.844 0.346 
1.900 2.988 0.345 
2.000 3.132 0.344 
2.100 3.276 0.343 
2,200 3.421 0.343 
2.300 3.566 0.342 
2.400 3.711 0.341 
2.500 3,857 0.341 
2.600 4.003 0.340 
2.700 4.150 0,340 
2.800 4.296 0.339 
2.900 4.443 0.339 
3.000 4.590 0.339 
3.100 4.737 0.338 
3.200 4.884 0.338 
3.300 5.032 0.338 
3.400 5.180 0.337 
3.500 5.328 0.337 
3.600 5.476 0.337 
3.700 5.624 0.337 
3.800 5.112 0.337 
3.900 5.920 0.336 
4.000 6.068 0.336 
4.100 6.217 0.336 
4.200 6.365 0,336 
4.300 6.514 0.336 
4.400 6.662 0.336 
4.500 6,811 0.336 
4.600 6.960 0.336 
4.700 7.108 0.336 
4.800 1.257 0.335 
4.900 7.406 0.335 


In der Tabelle sind optimale Werte der Windradberechnung gegeben. Die Werte wurden 


aus der Schnellaufzahl A und dem optimalen Abminderungsfaktor £opı = f(A) berechnet. 


Es ist: ® 


= X 
er kopt 


=c 


Czpl/r 


11.149 


56.193 
52.460 
48.867 
45.466 
42.290 
39.358 
36'672 
34.227 
32.009 


30.000 
28.183 
26.537 
25.046 
23.692 
22.460 
21.337 
20.310 
19.370 
18.506 


17.710 
16.976 
16.296 
15.666 
15.080 
14.534 
14.025 
13.549 
13.103 
12.684 


12.290 


&,pl/r = ozt/r 


1=% 


x 


5.000 
5.100 
5.200 
5.300 
5.400 
5.500 
5.600 
5.700 
5.800 
5.900 


6.000 
6.100 
6.200 
6.300 
6.400 
6.500 
6.600 
6.700 
6.800 
6.900 


7.000 


10.390 
10.539 


für e > 


0.589 
0.589 
0.589 
0.589 
0.589 
0.589 
0.590 
0.590 
0.590 
0.590 


0.590 
0.590 
0.590 
0.590 
0.590 
0.590 
0.590 


vuumubuuuvuu uw 


Czplir 


0.217 
0.029 
0.201 
0.194 
0.187 
0.180 
0.174 


Bilder 


ie] (umequrure?g) 
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Windkonverter - Bauarten 


Auftriebskonverter Widerstandskonverter 
horizontale vertikale horizontale vertikale 
Achse Achse Achse Achse 
Fe, | ? Kg Ez 
Langsam- gemäfigter Schnelläufer Darrieus-Rotor (freilaufender Savon.-Rot. diskrete kontinuierliche 
läufer Schnelläufer 1-, 2-, 3-blättrig (Urform ) Savonius-Rotor) im Widerst-Elem. Wid.- Elemente 
(X=2) (22 As#4) (A»4) induz. Wirbel 
Tandem- Mantel- tip-vane- alte Segel- Anemo- freilaufende Rotoren m. 
konverter turbine Konverter tuchmühle meter Rotoren Windleit- 






Sonderbauarten: 


| vorrichtung 


Savonius-Rotor Windturbine 


A.- Darrieus A. G.-Rotor 


a) primärwirbelinduzierte 
Ansaugkonverter H- Darrieus Herter-Rotor 





b) Heißluft- Auftriebskonverter Hansen-Rotor Savonius-Rotor Sav.-Rotor 
mit luvseitigem mit rad. 
Verdränger Leitblechkonz 






(Fl ügelelement) 


Bild 2.1: Radelement und Blattelement (Flügelelement) [10] 





Bild 2.2: Druckverlauf in der Nähe der Rotorebene 4) 





Bild 2.3: Strömung um ein Windrad [4] 





Bild 2.4: Strömungsgrößen an einem Flügelelement |15 
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Bild 2.5: Winkelbezeichnungen am Flügel [2] 


Konuswinkel 





Der Umlaufwinkel hat den Wert Null (8 = 0), wenn der Flügel seine 
tiefste Stellung eingenommen hat. Diese Stellung ist im Bild mit ”*” 
gekennzeichnet. # läuft - beim Blick mit der Anströmrichtung - in 


Richtung des mathematisch positiven Umlaufsinnes. 


Bild 2.6: Achswinkel r, Konuswinkel x und Umlaufwinkel 6 [15] 





Bild 2.7: Idealer Leistungsbeiwert c,; als Funktion der Schnellaufzahl A mit dem 


Abminderungsfaktor € als Parameter [3] 
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Bild 2.8: Idealer Leistungsbeiwert c,; als Funktion des Abminderungsfaktors £ 
mit der Schnellaufzahl A als Parameter [11] 
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Bild 2.9: Flügelwirkungsgrad np als Funktion der Schnellaufzahl A mit der Pro- 
filgleitzahl e als Parameter. (Für die Berechnung der Kurven wurde 
für jede Schnellaufzahl der optimale Abminderungsfaktor Eopı gewählt) 


a 





Bild 2.10: Flügelumrisse für das optimale Rotorblatt im Vergleich mit zwei Kom- 


promißformen [4] 
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Bild 2.11: Blattzahlverlustfaktor n, nach Hütter als Funktion der Schnellaufzahl 


A: mit der Flügelzahl = als Parameter [11] 
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Bild 3.1: Mit Panels belegte Oberfläche einer Flügel-Rumpf-Kombination [21] 





c 


Bild 3.2: Strömung um einen Tragflügel mit verschieden großer Zirkulation [25]: 
a) hinterer Staupunkt auf Oberseite 
b) hinterer Staupunkt auf Unterseite 
c) Staupunkt in Hinterkante, glatter Abfluß: Kuttasche Abflußbedingung 


erfüllt. 





gte ”Innere Fläche” einer Flügel-Rumpf-Kombination 


Bild 3.3: Mit Panels bele 





Bild 3.4: Darstellung des Gitternetzes zur Beschreibung der Körpergeometrie 
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Bild 4.1: Koordinatensystem am panelierten Flügel des Windrades Nr.7 der Ae- 


rodynamischen Versuchsanstalt zu Göttingen. 


Bild 4.2: Panelgrundverteilung in Richtung der Flügeltiefe nach Kap. 4.1 am 
Beispiel des Windrades Adler 25 (Profil am Radius r=7.86m). 
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Bild 4.5: Wirbelleiter, bei der der Bereich an der Flügelspitze dem Nabenbereich 


vorauseilt. 
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Bild 5.1: Druckverteilung an einem Profil NACA 2412; berechnet mit dem Pa- 


nelverfahren. 
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Bild 5.2: Gemessene Druckverteilung an einem Profil NACA 2412 [31]. 
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Bild 5.3: Druckverteilung an einem Profil NACA 0015; berechnet mit dem Pa- 


nelverfahren. 
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Bild 5.4: Gemessene Druckverteilung an einem Profil NACA 0015 [32]. 
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Bild 5.5: Profile am Flügel des Windrades Adler 25 [34]. 





Bild 5.6: Flügel des Windrades Adler 25, so wie er im Unterprogramm zum 
Panelverfahren modelliert wurde. Gezeigt ist nur die Panelverteilung in 
Flügellängsrichtung (15 Panels), wie sie für die Programmläufe gewählt 


wurde. 





Bild 5.7: Strak des panelierten Rotorblattes des Windrades Adler 25. 
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Bild 5.8: Übersichtszeichnung zum Windrad Nr. 7 der Aerodynamischen Ver- 


suchsanstalt zu Göttingen [35]. 





Bild 5.9: Strak des panelierten Rotorblattes des Windrades Nr.7 der Aerodyna- 


mischen Versuchsanstalt zu Göttingen. 
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Bild 5.10: Druckverteilung am Flügel des Windrades Adler 25 (Radius: r=9.53m). 
Die Druckverteilung wurde mit dem Panelverfahren für die 3dimensio- 
nale Strömung berechnet. X, = 6.2, &, = —1.2°, Ymaz = 90°, 6 Wir- 
belleitersprossen. 


a) 10 Panels pro Flügelelement 
b) 20 Panels pro Flügelelement 


c) 40 Panels pro Flügelelement 
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Bild 5.11: Druckverteilung am Flügel des Windrades Adler 25 (Radius: r=9.53m). 
Die Druckverteilung wurde mit dem Panelverfahren für die 3dimensio- 
nale Strömung berechnet. A, = 6.2, fı = —1.2°, Ymaz = 90°, 6 Wir- 
belleitersprossen. 

Bei 20 Panels pro Flügelelement wurden verschiedene Panelverteilun- 
gen in Flügellängrichtung gewählt: 

a) 10 Panels 

b) 20 Panels 





Bild 5.12: Gemessener Leistungsbeiwert cp des Windrades Adler 25 als Fuktion 
der Schnellaufzahl an der Flügelspitze A; mit dem Flügelwinkel £, als 
Parameter [34]. Für ausgewählte Punkte des Kennfeldes (DO) wurde 
der Leistungsbeiwert mit dem Panelverfahren nachgerechnet. 

Es bedeutet: 
V: Ymaz = 90°, 6 Wirbelleitersprossen, 
A: Ymaz = 360°, 24 Wirbelleitersprossen, 


®: Ymaz = 720°, 24 Wirbelleitersprossen. 
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Bild 5.13: Vergleich von 3dimensionaler und 2dimensionaler Rechnung mit dem 
Panelverfahren 
am Profil FX 77 W 343 (Radius r=3.75m) des Windrades Adler 25. 
a) Druckverteilung aus 3dimensionaler Strömungsberechnung am Wind- 
rad. A; = 6.2, A = —1.2°, Ymaz = 90°, 6 Wirbelleitersprossen. 
b) Druckverteilung aus 2dimensionaler Strömungsberechnung am Pro- 
fil bei einem den Windradverhältnissen entsprechenden Anströmwinkel 
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Bild 5.14: Vergleich von 3dimensionaler und 2dimensionaler Rechnung mit dem 
Panelverfahren 
am Profil FX 77 W 258 (Radius r=5.45m) des Windrades Adler 25. 
a) Druckverteilung aus 3dimensionaler Strömungsberechnung am Wind- 
rad. Aı = 6.2, ßı = —1.2°, Ymaz = 90°, 6 Wirbelleitersprossen. 
b) Druckverteilung aus 2dimensionaler Strömungsberechnung am Pro- 
fil bei einem den Windradverhältnissen entsprechenden Anströmwinkel 
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Cp(-) 








Bild 5.15: Vergleich von 3dimensionaler und 2dimensionaler Rechnung mit dem 
Panelverfahren 
am Profil NACA 4421 (Radius r=6.85m) des Windrades Adler 25. 
a) Druckverteilung aus 3dimensionaler Strömungsberechnung am Wind- 
rad. A; = 6.2, 8, = —1.2°, Ymaz = 90°, 6 Wirbelleitersprossen. 
b) Druckverteilung aus 2dimensionaler Strömungsberechnung am Pro- 
fil bei einem den Windradverhältnissen entsprechenden Anströmwinkel 
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Bild 5.16: Vergleich von 3dimensionaler und 2dimensionaler Rechnung mit dem 
Panelverfahren 
am Profil NACA 4418 (Radius r=8.75m) des Windrades Adler 25. 
a) Druckverteilung aus 3dimensionaler Strömungsberechnung am Wind- 
rad. A: = 6.2, fı = —1.2°, Ymar = 90°, 6 Wirbelleitersprossen. 
b) Druckverteilung aus 2dimensionaler Strömungsberechnung amı Pro- 
fil bei einem den Windradverhältnissen entsprechenden Anströmwinkel 


a = 4.6°. 
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Bild 5.17: Vergleich von 3dimensionaler und 2dimensionaler Rechnung mit dem 
Panelverfahren 
am Profil NACA 4415 (Radius r=12.25m) des Windrades Adler 25. 
a) Druckverteilung aus 3dimensionaler Strömungsberechnung am Wind- 
rad. A, = 6.2, 3, = -1.2°, Ymaz = 90°, 6 Wirbelleitersprossen. 
b) Druckverteilung aus 2dimensionaler Strömungsberechnung am Pro- 
fil bei einem den Windradverhältnissen entsprechenden Anströmwinkel 
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Bild 5.18: Auftriebsverteilung über dem Windradflügel. Aufgetragen ist der Auf- 
triebsbeiwert c„ über dem normierten Windradradius. Die stark aus- 


gezogene senkrechte Linie kennzeichnet den normierten Radius der 


Flügelwurzel. 





Bild 5.19: Gemessene Windradkennwerte (Leistungsbeiwert cp, Momentenbeiwert 
cm und Widerstandsbeiwert cw) des Windrades Nr. 7 der Aerodyna- 
mischen Versuchsanstalt zu Göttingen als Funktion der Schnellaufzahl 


an der Flügelspitze A, [36]. Die mit dem Panelverfahren berechneten 


Werte sind mit ”x” eingezeichnet. 
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Bild 6.1: Windradkennzahlen: idealer Leistungsbeiwert cp;, Flügelwirkungsgrad 


np, Blattzahlverlustfaktor n, und Gesamtleistungsbeiwert cp berech- 
net nach Hütter als Funktion der Schnellunfzahl an der Flügelspitze A, 


für das Windrad Adler 25 bei einem Blattwinkel 8, = —1.2° 
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Bild 6.2: Gesamtleistungsbeiwert cp als Funktion der Schnellaufzahl an der Flügel- 


spitze A, mit dem Blattwinkel 8, als Parameter berechnet nach Mütter 


für das Windrad Adler 25. 
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Bild 6.3: Gesamtleistungsbeiwert cp als Funktion der Schnellaufzahl an der Flügel- 


spitze A, bei verschiedenen Flügeltiefen. 
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Bild 6.4: Gesamtleistungsbeiwert cp als Funktion der Schnellaufzahl an der Flügel- 


spitze A, bei verschiedenen Flügelzahlen aber konstanter Größe des 


Wertes z-t. 


Anhang 


Herleitung von Gl. 2.12 


Es soll ein Flügelement dA =t- dr betrachtet werden. Auftriebs- und Widerstandskraft 


an einem Flügelelement sind 


und 


w ist die Relativgeschwindigkeit am Flügel, die Geschwindigkeit, der der Flügel ausgesetzt 
ist: 


w=./vU+u : (1) 
In Axialrichtung wirken die Schubkomponenten 
Fıs = Fa:cosp 
und 
Fw,s = Fw - sinp 


Aus einer Betrachtung eines Radelementes bei einem Windrad mit x Flügeln erhält man 
als Gesamtschubkraft 
Fs=z(Fas+Fws) - (2) 


Als nächstes wird die Gesamtschubkraft auf das Radelement aus einer Betrachtung der 


Strömung berechnet. Es ist 


Fs = m(vı - v5) (3) 


m = p-?2rr:-dr:.vz(vı — v5) 


Setzen wir nun Gl.2 mit Gl.3 zusammen, so erhalten wir 


w’t(ca c0oSp + c„ sinp) = ?rrv (vı — v5) 


D|n 


Beachtet man, daß 


u2 
cosyp = — 
w 
und 
2 v2 
sing = — 
w 


ist, und nimmt man Gl. 1. Def. 2.2, Def. 2.4. Def. 2.5 und Def. 2.6 zu Hilfe, so erhält 


man nach einigen Umformungen Gl. 2.12. 
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